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THEORY OF GASTRIC FUNCTION 


J. TAMARIT TORRES AND F. ENRIQUEZ DE SALAMANCA 
THE UNIVERSITY OF MADRID 


On the basis of some assumptions concerning gastric emptying a 
theory of gastric function is developed which explains the known facts 
—chiefly the emptying of liquids according to an exponential law. On 
the basis of this theory an equation is derived which relates the rate of 
secretion to the quantity of liquid retained in the stomach after a cer- 
tain time. A study is made of the general characteristics of this equa- 
tion. Some assumptions concerning different formulae describing the rate 
of secretion are found useful in experimental investigations, particular- 
ly in Hollander’s dilution indicator method. 


One of the authors (Salamanca, 1943) has developed a method 
which allows us to determine the amounts of tea and gastric juice 
retained at any time after the ingestion of a quantity of tea. The 
continued use of this procedure over a period of many years has re- 
sulted in a great number of observations which require a theory of 
the gastric function for their correct interpretation. Such a theory 
is developed in the present paper. 


Schematically the function of the stomach consists in merely re- 
ceiving a certain alimental volume (in exploratory tests, 250 cc. of 
tea) and in emptying it slowly, mixed with a gastric secretion, which 
contains hydrochloric acid, chlorides and ferments as essential ele- 
ments. In practice, in order to characterize the composition of gas- 
tric juice, it is sufficient to determine the concentration of hydrochloric 
acid and of total chlorine; the determination of the ferments in the 
contents is not usually made in the clinic. 


We denote the volume of tea retained by the stomach by a func- 
tion of time, R(t), the quantity of gastric juice secreted between the 
time of ingestion and time ¢ by f,(t) and the juice retained by the 
stomach in time ¢ by f.(¢). 

We will first study gastric emptying, which is the result of stom- 
ach motility, and afterward the secreting function, distinguishing be- 
tween total secretion and secretion of the specific elements, hydro- 
chloric acid and total chlorine. 
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I. Theory of Gastric Emptying. 
The two initial hypotheses of our theory are the following: 


a) The pylorus receives a stimulus every wu seconds to which it 
is able to respond in two ways—opening or not opening. The first 
response has a probability p and the second the probability q, 
(po 1). 

b) Each time the pylorus opens, the stomach empties a constant 
fraction c of its contents. 


With these two hypotheses we are able to characterize the func- 
tion R(t). Let us suppose that 250 cc. of tea have been taken in. After 
the first opening of the pylorus, there will remain in the stomach 


250 — 250 c= 250(1—c) cm’ of tea. 
After the second opening there will remain 
250 (1 — c) — 250(1 — c)c = 250 (1— c)?2, 


and so on. After the nth pyloric opening there will be left in the stom- 
ach a quantity of tea 


R(t) =250(1—c)*. (1) 


When 7 becomes sufficiently large with respect to c, (1) will be ap- 
proximately equal to 


R(t) = 250 e-, ; (2) 
and as ” is proportional to t we may finally write 
R(t) = 250 e**, (3) 


Thus R(t) is a decaying exponential function of time. This has been 
confirmed by experiment in extractions made within 30, 60, 90, 120, 
150 and 180 minutes after the intake of a test meal (Salamanc¢a, Gar- 
cia Morato and Delgado, 1945). The value of the constant k is 0.0284 
for a normal adult and 0.0218 for a normal child (Salamanca and 
Tamarit, 1948a). 

Hypothesis a) has also been verified by experiment. The pyloric 
opening was observed radiologically (Tamarit and Truchuelo, 1948). 
Using a reasoning similar to the one used by E. Marsden and T. Bar- 
rat (1911) in the study of the emission of corpuscles by radioactive 
atoms, it is possible to calculate the probability that the interval be- 
tween two consecutive pyloric openings will be greater than a given 
interval 7’. Actually the pylorus will receive during this time a num- 
ber of stimuli equal to e(7) and the probability sought is approxi- 
mately 
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P= e-?Pé&(t) 
> 


where p , it will be recalled, is the probability that the pylorus opens 
upon receiving a stimulus. 

In the paper of J. T. Tamarit and J. Truchuelo (1948) an inter- 
esting fact was encountered, namely, the existence of a refractory 
time (p) in pyloric functioning. Thus after each pyloric opening there 
is a time of p seconds during which no stimulus reaching the pylorus 
can be effective. Taking the existence of this refractory phase into 
account, the number of stimuli reaching the pylorus during time 7 
will be 


op) / Bs, 
and the probability sought is 
P. = e?\t-p)/u . 


The theoretically calculated probability agrees satisfactorily with 
the values actually found. In the above observations it was found that 
for the normal stomach the following values hold 


p/u=5.91 X 107 sec; p= 5.14 sec. 


The average number of pyloric openings (n) occurring in time 
t is 


— Z it eeOG. 
u/p— p 

Thus if the tea meal is tested in different persons and the extrac- 
tion of the gastric contents is done at the same interval of time after 
ingestion, the quantities of tea retained will depend chiefly on the 
number of times the pylorus has opened. The quantity of the tea re- 
tained must therefore be related to the distribution of frequencies of 
the number of favorable cases of a phenomenon (pyloric opening) 
having a small probability (p), ie. according to Bortkiewicz-Pois- 
son’s law. This has been confirmed by J. T. Tamarit (1947) through 
the results of investigations carried out on adults and children. 

As to the value c , which might be called “emptying gastric pow- 
er,” there is to date no definite method for determining it. Assum- 
ing the rate of pyloric opening to be the same for tea as for barium 
substance used for roentgenologic observation, the value of ¢ might 
be calculated indirectly. One thus obtains the value c= 0.0104, that 
is, the stomach empties 1.04 per cent of its contents upon each caloric 
opening. Further investigations have been undertaken. 


7 
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II. Theory of Gastric Secretion. 


The first problem is to find the relation between quantities se- 
creted and quantities retained by the stomach. 

We assume that gastric emptying is performed according to 
equation (3). Differentiating equation (3) we have 


dR (t) 
dt 


— k- 250 e* =—k-R(#). 


Now we shall consider the volume of secreted liquid. 


Ill. Secretion of Gastric Juice. 


According to our notation the quantity of gastric juice retained 
in the stomach during the time t is f,(¢). In an infinitesimal interval 
of time, dt, the variation in the amount of gastric juice retained, 
which is f,'(t) dt, will be equal to the volume secreted in this inter- 
val [f,'(t)dt] less the quantity of juice that has been emptied 
[kK - fo(t) - dt]. We may thus postulate the following differential 
equation: 


fe’ (t) =f,’ (t) — kf (¢). (4) 
This gives 


f(t) =e" | oF f(t) dt + c | 5; (5) 


where C is a constant of integration. Imposing initial conditions so 
that f2(¢) = Jo fort = 0 (Jo = gastric juice contained in the fasting 
stomach), we have 


f(t) =e" | {f ek f(t) dt + J, (6) 


0 


The solution of the integral of (6) depends on the form of the 
function f, (t). 


In the simplest case f,'(t) = constant = VY», and elueun (6) 
takes the form of 


(Xt) = voCL—e*") 7k Seer, rye 


In this case the quantity of juice retained in the stomach will 
increase with time approaching a limit, 


Lim. f.(t) =v,/k. | ape (8) 


t=00 
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Thus if the stomach secretes at a constant rate, an equilibrium 
condition will be reached with the quantity of retained juice equal to 
Vo/k. At this time the secretion rate and the emptying rate are equai. 

This is what occurs in the resting stomach during fasting. Ac- 
cording to K. Faber (1935), one can extract from a normal person a 
quantity of juice of approximately 20 to 50 cc., which corresponds 
to a resting secretion rate of .57 to 1.42 cc. per min. respectively. We 
have calculated an average value of 1.07 cm*/min. for adults and 
0.21 cm*/min. for children (Salamanca and Tamarit, 1948b). 

If the condition of gastric rest is altered by a stimulus (tea or 
test meal, histamine, insulin, etc.), we have a different situation. 
However, even under such circumstances one may proceed within 
short intervals of time as if the rate of secretion were constant, as 
some authors have done (Hollander and Penner, 1940; Schoen and 
Knoefeld, 1947; Schoen and Griswold, 1947). Their results ought to 
be considered only as first approximations, especially because they 
make the additional hypothesis that the rate of emptying is also con- 
stant during the period of observation. With our theory one obtains 
a better approximation to reality. 

Before making any complementary hypothesis on the form of the 
function f,'(¢) in the active stomach, it is most convenient to study 
the conditions of maximum and minimum of the function f,(¢). Dif- 
ferentiating equation (6) we have: : 


f.' (t) =— kf. (t) + fi’ (t) ; (9) 
fo.” (t) =— kf,’ (t) + f," (€) =F? f.(t) — kf,’ (t) + fi" (t). (10) 
In order that f,(t) has a maximum, we must have 


f(t) = kf. (t) 
and also 
fo” (t) =f," (t) <0. 
Similarly for the amount of juice retained to be a minimum, we must 
have f,'(t) = kf2(t); fe" (t) =f," (t) > 0. 

Experimentally it is not possible to make direct determinations of 
f,/(t) and we only have a direct knowledge of f.(¢) ; thus the exper1- 
mental observations are useful to deduce consequences about f,’(t) 
from the known values of f.(t). This function has been studied ex- 
perimentally by F. Enriquez de Salamanca and J. Picazo (1948) in 
adults and by F. Enriquez de Salamanca, V. Garcia Morato and C. 
Perez Delgado (1945) in children. From these observations one de- 
duces that the function f.(¢) at first increases, attains its maximum 
within 60 minutes, then decreases to attain a minimum within about 
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120 minutes, and afterward it increases again. No observations have 
been made in times greater than 150 min., but it is plausible to suppose 
that the function will have another maximum and afterward will di- 
minish asymptotically to again attain its resting value. In children 
the function f,(t) behaves similarly with the difference that the mini- 
mum is attained within 150 min. The maximum value of the juice re- 
tained in adults is about 168.33 and in children 28.57 cc. 

According to the above conditions, necessary for the existence of 
turning points of the function f.(t), we have that within 60 min. the 
rate of secretion has to be 4.78 cc. per min. (= 0.0284 X 168.23) and 
in decreasing phase. Since the resting rate is 1.07 cc. per min. it must, 
of course, have passed through a maximum before the 60 min. With- 
in 120 min. in adults the secretion rate must be 2.54 cc. per min. 
(= 0.0284 X 89) and in increasing phase; thus it must have passed 
through a relative minimum. After the 120 min. the secretion rate 
will continue increasing, it will pass through another maximum and 
then it will diminish tending asymptotically toward the value at rest. 
Similar deductions can be made referring to children. 

In short, the conditions imposed by observations on the function 
f,'(t) during the digestive process are: 


1. The function f,’(¢) tends to approach the value of the secre- 
tion rate at rest, 


Lim. f1'(t) =. 


t=0o 


2. The function f,’(¢) must be such that it makes f.(t) twice 
maximum once within 60 min. and the second time within 150 min. 
more or less. It must also have the property to make f.(¢) minimum 
at about 120 min. 


3. The function f,'(¢) must have a maximum (absolute or rela- 
tive) between 0 and 60 min., a relative minimum between 60 and 120 
min. and a relative maximum at about 150 min. 


We believed this digression necessary first as an orientation for 
the empiric investigator (physiologist, experimental pathologist) and 
secondly as a point of departure for new researches that might be 
undertaken to resolve important questions of the gastric physiology. 

In the following some interesting hypotheses are considered in 
relation to the imposed conditions. 


First Hypothesis. Let us suppose that 
fi’ (€) = Vo + vie! (11) 
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and, in the most common case, let us suppose that k, be distinct from 

the constant of emptying, which we will designate by k,. The gastric 

response consists actually in these cases of a sudden increase followed 

by a slow decrease, and discounting the first short phase, the rest may 

be represented by an exponential function of the type (11). A simi- 

lar thing is done by N. Rashevsky (1938) studying the delayed reflex. 
It can easily be seen that this hypothesis satisfies condition 1. 


Substituting (11) into (6) we have 


fo(t) = Uo/ky + v1 (et — et) / (ke, — ke), (12) 
and its derivative is 
fo’ (t) = v, (h,e™! — ke") / (ky — ke), (13) 


which vanishes for the following value of the time 


a log k, — log kz 
(ky a Ke) 
Since f,”(¢) < 0, in that moment there is a maximum for f.(t) and 
thus condition 2 is partly satisfied. In order that the maximum coin- 
cides with 60 min, we must have, taking into account the value of k, 
in normal adults, k, = 0.0083. 

With this first hypothesis the moment of the maximum of f2(t) 
does not depend upon v,, but on the constants k, and k, according to 
formula (14). The time 6 is equal to 60 min. when k, = 0.0083; it 
decreases with increasing k, and has a minimum when k, = k,, then 
it increases again. 

Substituting (14) into (12) we find that in the maximum the 
quantity of gastric juice retained will be 


f2(0) =Vo/hy + (vie™") /Ke ; (15) 


a value that depends directly on v,. In order that in a normal adult 
the maximum of juice be 170 cc. of gastric juice retained within 60 
min. it is necessary that v, be equal to 6 cm?/sec. This would be the 
initial secretion rate according to (11) and its highest value in the 


course of the digestive process. ! 
The case in which k&, and i» are equal deserves special considera- 


tion. It is easy to see that 
fo(t) =v (1+ €*)/k + vte™ , (16) 
and the maximum of f.(t) is now attained in the time 
b= 1/k. 


(14) 
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This value 6 is also obtained putting k, = k, into formula (14) 
and applying l’Hopital’s rule. In this case the maximum would be 
attained, taking the normal value of k, into account, within 35.21 min. 
in the normal adult and within 46.95 in the child. 

This hypothesis does not satisfy all of the above conditions; it 
fails because it does not make f.(t) minimum within 120 min. nor 1s 
there-a maximum after 150 min. 


Second Hypothesis. Let us suppose that 
(t-t1)2 (t-te)? 


Fy, 2012 2 2022 
f,' (t) =v. + ——e +—e ; (17) 
V 2a V 2a 


According to this hypothesis, the secretion rate is represented 
by two waves superposed upon the resting secretion rate, each of 
which is a normal Gauss curve. This hypothesis is plausible on the 
following grounds. The secretion of gastric juice is the result of the 
activity of several glandular units. It may be supposed that they do 
not all function simultaneously, but that the number acting in any 
interval of time is expressed by the normal probability function, and 
the quantity of juice secreted is proportional to the number of units 
acting. The factors F, and F, in this hypothesis are related to the 
total number of secretory units present in the stomach and to the se- 
cretion rate of each isolated secretory unit. 

In our analysis we assume for simplicity that o, = o.. In what 
follows we show that the conditions imposed upon the function f/f,’ (¢) 
are satisfied with this hypothesis and the facts established by the em- 
piric observation are accounted for. 


Putting (17) into (6) we have 


F-. B1 e-kt 2 F. B2 e-kt 2 
A f ere 2 ff evrde, (18) 
V2n 20 V 2a Zo’ 


fo(t) =Vo/k + 


where 
2= (t—a,)/o; a —t,'+ ko; b= k(t,'+ a) /2; 


Z= (€— a2) /o; Ot, = ts + k o?; fo = k(t + 2) /2. 


It is easy to see in (18) that if ¢ increases indefinitely the last 
terms of the sum tend to vanish, so that for t = 0, f2(t) = v./k, 
i.e., condition 1 is fulfilled. 

As to condition 2 it is always satisfied when t, and t are suf- 
ficiently separated so that in the first moments the secretion rate is 
only expressed by the first two terms of (17) and afterward by the 


(19) 
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first and third. Referring thus to the beginning of the digestive pro- 


cess we have 
eet Se 
eal e-**/2 dz — =a) e-*?/2 ic| , (20) 
V 2a -00 V 2n ae 


where we may neglect the second integral as z is always greater than 
2. Moreover if we replace the first integral by the tangent at its in- 
flection point we have 


fo(t) =v,/k + F, e* e& [0.3989 (t —a,) + 0.50]. (21) 


Making the derivative of f.(t) equal to 0 we find that this func- 
tion attains its maximum in the moment 


6=1/k + a, — 1.253 o (22) 


f.(t) =v/k+ Roches | 


if the maximum coincides with 60 min. Then we have 
ko? —1.2530+ t,'+1/k—60=0. (23) 
This equation gives us the relation that ought to exist between 


the values ft, and «. For each value of t, two values of o are solu- 
tions of (23). In order that the values of o be real it is necessary that 


4.3782 —4kt,20, 


i.e., that t, be less than or equal to 38.62 min. We shall only consider 
the case in which t, = 38.62 which gives one value foro. In that case 
we find o = 22 min. Substituting these values into (19) we have 


a, = 52.45; 6, = 1.298. 


Replacing the two known values in (18), putting ¢ = 60, at 
which time f.(¢) has its maximum, ignoring the second integral, and 
calculating the first between z) = —2.38 and z = 0.32 by means of 
A. L. Bowley’s tables (1920) we obtain 


170 = 87.7 + 9.35 F3; Fi,=14.1. 


We also may calculate F, if we substitute the known values in 
(17) bearing in mind that f,'(60) = 4.83. 
4.83 = 1.07 + 0.2492 F,; Fi 1b.) . 
The difference between the values of F, is due to our approxi- 
mation. We shall take for the value of F, the average of the two re- 
sults = 14.6. 
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To continue our analysis in the vicinity of 120 min. we may de- 
velop f.(¢) in Taylor’s series, bearing in mind that within 120 min. 
the function f,(t) has a relative minimum (= 89 cc.) and thus, as 
shown before, f2'(120) = 0 and f.” (120) =f,” (120). Including in the 
development only the terms of the second degree, we have 


f2(€) = f2(120) + fa" (120) - (6— 120) 2/2; (24) 


ie., in the vicinity of 120 min. the function f,(t) may be represented 
by a parabola. Keeping in mind that within 60 min. the value of f. (t) 
is 170 cc. we are able to calculate 


f.” (120) = 0.045. (25) 


Putting (24) into (6), introducing a slight change and differen- 
tiating the two members, we find a fairly exact expression of Tanke ye 
valid in the vicinity of 120 min. 


eEfs(120) + f," (120) - (t—120)2/2] = fi ot f(t) dt + Jo: 
ke* [f, (120) + f," (120) - (t—120)2/2] (26) 
+ o f," (120) - (t —120) =e f,'(t) : 
fi (t) = kf (120) '+ f.” (120) - (— 120) 
+ kf" (120) - (t—120)2/2. 


This expression of f,'(t) tells us that the secretion has a mini- 
mum in the time 


7 = 120 —1/k = 84.79. 


The value of ¢, may be calculated in the following way: We know 
that f,”(¢) at 120 min. has a value of about 0.045, calculated in (25) ; 
the secretion rate in that time is 2.54 and the secretion rate, due to 
the first wave calculated by Charlier’s tables is 0.02. Thus the se- 
cretion rate due to the second wave is 1.45 ec. (1.07 + 0.02 + 1.45 
= 2.54). If we differentiate (17 ) and denote the exponentials by Vi 
and V. respectively, we have 


fi" (¢) =— (€—t) - Vi/o? — (t— t,) - Vo/o?. (27) 


Substituting the known values when t = 120 we have: 0.045 = 
— 1.627/484 — 1.54(120 — t,)/484, which gives t, = 135. There- 
fore in order that V. be equal to 145 within 120 min. it is necessary 
that F, = 4.58. Thus the function (17) is completely characterized. 


After 120 min. in the first integral of (18) z is greater than 3 
and we have 
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fo(t) =v./k + F, 8: e*t g 
; Ee (28) 
+ FF, e8 e**[0.3989(t—a,) + 0.56], 
which is a good approximation for times greater than 120 min. It is 
easy to see now that after 120 min. the function fo(t) will pass 
through a second maximum in the time 
ve o ef 
6= ao, + 1/k — 1.253 « — —_—_—_—_—_ 29 
: ” 0.3989 F, ef ee 

which, when the known values are substituted, gives 145 min. 

The value attained by f.(¢) in this second maximum may be 
calculated by substituting the value of (29) in (28) 


fo(0) =V./k + 0.8989 F, ef e-*°/k . (30) 


Putting the known values of the constants into this equation, we 
find that the quantity of juice retained in the stomach at this second 
maximum must be approximately equal to 97 cc. 

This completes the analysis of our second hypothesis. It explains 
the known facts and fulfills the conditions imposed on the function 
f,'(t) . There remains the verification of the relation between t, and 
o, expressed by (23). 


Third Hypothesis. In this hypothesis we unite the two former 
ones, i.e., we suppose that 
(t-t1)? 
fi(t) =Vo +, 6 + (Fy/V2aje 2 
(t-te)? 


ie (Fo/V2a)e 2 


In this hypothesis the gastric juice secreted during the digestive 
process is the result of the mixture of four different secretions: (a) 
resting secretion (Vv); (b) secretory fraction of exponential type; 
(c) two secretions superposed on the two former ones in the form of 


two Gauss’ curves. 
Putting (31) into (6) we have 


fi (t) = 0./h, Fv, (oe — e*) / (es — tee) 


z z? 
Bs eat 
of EEE are ad a [ e*/2 dz, 
V2n ; Van oh, 


20 


(31) 


(32) 


with the same substitutions as in (19). Referring to the beginning 
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of the digestive process we may neglect the second integral and sim- 
plify the first, as in the former hypothesis. Thus we have 


fo (t) = Vo/K'+ V1 (e™* — e™*) / (les — Ia) 


(33) 
+ F, et e&[0.8989(t—a,) + 0.50]. 
Putting its derivatives equal to zero gives 
f.' (t) =v (ka es! — Ke eat) / (ky — k,) + 0.8989 Br, et eb: ( ) 
34 


— k, F, e& e'[0.3989 (t—a,) + 0.50] —0. 


Dividing both members by e*' and taking only the first two 
terms of the following development 


eVr-ks)t — 1 + (kh, — hn) tes--, (35) 
we get 
V,— V1 ke t — 0.8989 F, k, c(t —a, + 1.253 0—1/k,) =0. (86) 


This relation is verified only in the moment of the maximum, i.e., 
when t= 60. The quantity a, — 1.253 « + 1/k corresponds to the time 
in which the maximum of f.(t) is attained in the former hypothesis. 
We shall denote it by 7 and denote the factor 0.3989 F,k,e® by h. We 
finally obtain 


v,=h(60—T)/(1— 60 k). (37) 


From this equality one may deduce a consequence for k.. If 
v, > 0, and if T < 60, we must have k, < 0.0117; if T > 60, then 
e001 

As in the case of the former hypothesis, we now have the follow- 
ing equation to study the relations between ¢, and o 


Kyo? — 1. 2b3.00F 1/i + & TO . (38) 
The condition for the vanishing of the discriminant is 
T —t,= 21.5 (39) 


in which case o = 22. The value of t, remains undetermined. Since 
in practice one admits the normal curve to have the characteristic of 
comprising nearly the totality of the observations between — 2 and 
+ 2o, we assume, in order to take the earlier activity of the glands 
during the digestive process into account, that t; = 2c. Any other 
hynothesis may be made a priori giving rise to new orientations for 


emmric investigations that may justify with a greater exaciivude the 
relation between t, and o. 
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Thus we have obtained the following values: 
t, = 44; T= 65" a, = 57.75; B1= 1.445. 


To calculate an approximate value of k, we shall use the same 
device as in the former hypothesis. We shall now replace the func- 
tion f2(t) by a parabola, 


f(t) = f2(60) + f,” (60) - (¢—60)?/2. (40) 


Remembering that for t= 0, f.(t) = 37.7 cc., we may calculate the 
value f,” (60) =— 0.0735. 
Putting (40) into (6) and proceeding as before we find 


fi (t) = ky fo (60) + f,” (60) 
- (t—60) +k, f.” (60) - (t— 60)?/2, 
which shows that the function f,'(t) passes through a single maxi- 
mum point at the time 7 = 60 — 1/k, = 24.79. 
That is to say, the function f,’(t) has a maximum at about 
25 min., therefore f,” (t) must vanish at this moment. Thus differen- 


tiating (31), ignoring the last term of the sum, and setting the deriva- 
tive equal to 0, we get 


(41) 


(ae (Ra Fy 202 


ie e =0; (42) 
o 7 


Tia (7) =k, Vy OE 


Substituting the known values and replacing v, by its value giv- 
en in (37) we have 


0.0108 F, = 0.3989 F,, e® k, +5 hk e?/ (60 k, —1), (43) 
and substituting approximately 
e252 —= 1 — 25 ky----, 
we finally obtain 
3.8 k.? + 0.4968 k, — 0.0108 =0, (44) 


which gives k, = 0.019. The other root of k, of the equation (44) is 
negative and must be discarded, since otherwise we would have a se- 
cretion which would increase indefinitely in the course of the diges- 


tive process. 
Substituting the value of k, in (37) we have 


V,= 1.082 F,, 


and putting the same value into (42) results in 
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inco—k OOM G59) he, 


We obtain two slightly different coefficients because of the ap- 
proximate character of some of the substitutions used, but the differ- 
ence is so small that the reasoning followed is quite justified. We may 
take the average value, which is approximately unity, i.e., 


Vg ee 


The values of v, and F, can now easily be calculated since they 
have to be such as to give f.(t) a value equal to 170 cc. within 60 min. 
and f,(t) a value of 4.83 cc./min. in the same time. Substituting the 
known values in (31) and (32) we have 


3.76 = 0.6313 F, , from which F; = v,— 5.91; 
132.30 = 23.511 F,, , from which F, = v, = 5.63. 


Taking the average value, we have Ff, = v, = 5.77. 

Now we can calculate each of the values of the fraction of secre- 
tion at 120 min. and we find that the value due to the second wave 
is 0.87 cc./min. By the same process as in the former hypothesis we 
now calculate t, which turns out to be equal to 152 min. The quantity 
F, = 6.23 , slightly greater than F, , but is approximately of the same 
order. 

Herewith we complete the analysis of our third hypothesis, which 
also explains the known facts of the gastric physiology and which 
fulfills the conditions imposed on function f,(t). The analysis of hy- 
drochloric acid and chlorides will be made in a later paper. 

This work has been carried out under the auspices of the Insti- 
tute for Experimental Medicine of the Higher Council for Scientific 
Research in Spain. 
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A society composed of individuals each of whom can perform one of 
two mutually exclusive activities R, and R, is considered. The tendency 
toward the performance of those activities is measured by the inten- 
sities ©, and «, of excitation of two corresponding neural centers, which 
cross-inhibit each other. It follows from the theory developed by H. D. 
Landahl that an individual with «, — *, — 0, that is one who has no 
preference for either one of the two activities, will on the average per- 
form FR, and R, with equal probability. As €, — &, increases, the prob- 
ability P, of R, increases, tending to 1. As «, — ®, increases, the prob- 
ability P, of R, increases, tending to 1. We have P, + P, == 1. 

The effect of imitation is now studied. The total number of individ- 
uals in the society which exhibits an activity R, at a given time is con- 
sidered as constituting a stimulus which increases ¢,. Similarly, the to- 
tal number of individuals which exhibits activity R, at a given time con- 
stitutes a stimulus which increases «,. Using the standard equations of 
the mathematical biophysics of the central nervous system, equations 
are established which govern the behavior of such a society and the fol- 
lowing conclusions are reached. 

If the quantity «, — «, is distributed in the society in such a way 
that the distribution function is symmetric with respect to «, — *, — 0, 
then on the average one-half of the population exhibits R,, the other 
half R,. This social configuration may, however, be unstable. The slight- 
est accidental excess of individuals exhibiting, for example, R,, may 
bring it into a stable configuration, in which most individuals exhibit 
R,, and only a smaller fraction exhibit R,. A slight initial deviation 
in favor of R, brings it into a stable configuration, in which most indi- 
viduals exhibit R,. Thus in this case there may be two stable configura- 
tions. If the population is in one of those stable configurations, and the 
distribution function of «, — #, is made asymmetric, favoring the other 
activity, the population will pass into a stable configuration, in which 
that other activity is predominant, if the asymmetry of the distribution 
exceeds a threshold value. ; 

By making some drastic simplifications the equations derived here 
may be reduced to a form which was postulated by the author previously 
in his mathematical theory of human relations. 


Consider a population of N. individuals. Let each of them have 
a tendency to perform two mutually exclusive reactions R, and R2. 
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These tendencies are measured by the corresponding excitations ¢, 
and «, which lead to R, and R., and which also cross-inhibit each 
other, as in the neural circuits introduced originally by H. D. Lan- 
dahl (1938) and used in our preceding paper (Rashevsky, 1949). 
Whether the reaction R, or R. is made is determined by the quantity 
€, — & (Landahl, 1938; Rashevsky, 1948b). Therefore we shall con- 
sider this quantity ¢ = e, — «,. An individual with ¢ = 0 has an 
equal tendency to R, and R, , and the probability of either one of them 
is 4. For ¢ > 0, the individual has a preference to R,. The prob- 
bability P. of R, tends to zero as ¢ tends to +a. For ¢ < 0 the 
preference is in favor of R,, and the probability P, tends to 1 as ¢ 
tends to —o . Let the number of individuals whose ¢ lies between ¢ 
and ¢ + d¢ be given by N(¢)d¢. We have 


[-N@)ag=N. (1) 


If we consider a sufficiently large number of individuals all with 
the same ¢, then at any time a fraction P,(¢) of them will perform 
reaction R,, and a fraction P;(¢) = 1 — P,(¢) will perform reac- 
tion R,. Hence, denoting by x(¢)d¢ the number of individuals with 
a ¢ in the interval between ¢ and ¢ + dd, and performing reaction 
FR, , we have 


a(¢) =P;($)N(¢). (2) 


If y(¢)d¢ denotes the number of individuals in the interval ¢, 
¢ + dd which performs reaction R, , we have 


y (6) =P2(¢)N(¢). (3) 


Now we shall consider the effect of imitation of behavior. When 
an individual imitates others, this may be interpreted in the following 
way: 
Each other individual performing R, acts on the given individual 
as a stimulus which adds to the excitation «,. Each other individual 
who performs reaction FR, acts on the given individual as a stimulus 
which adds to e.. The additional amounts ¢,' and 2’ vary according 
to the usual differential equation (Rashevsky, 1948b) 


OSes 4 
an ae. (4) 
We shall consider here the simple case in which the constants A and 
a are the same for all individuals, regardless of their ¢. The more 
general case leads to great mathematical difficulties. 
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Denoting by X the total number of individuals who perform R,, 
and by Y the total number who perform R., we have 


de,’ 
spe a eg (5) 
ee ay et (6) 
dt | 
Putting 
& —& —y, (7) 
we have 
aa yay (8) 
dt 
If the distribution function N(¢) is symmetric with respect to 
@ = 0, then, in the absence of external influences, X = Y, or 


X— Y=0. Half of the individuals will perform R, , the other half 
k.. Therefore, according to equation (8), y will remain zero. 

However let an accidental fluctuation result in a slight excess of 
X. Then y will begin to grow. The quantity which now determines 
whether FR, or R, occurs is no longer $¢, but ¢ + yw. As ¢ + y in- 
creases, the probability P, for each individual to perform R, increases. 
Therefore X increases while Y decreases, so that X — Y increases. 
Similar considerations hold for an accidental excess of Y. Thus the 
configuration X — Y for a symmetric N(¢) may in general be un- 
stable. 

To investigate whether a stable configuration exists and to de- 
termine it, we must make specific assumptions about N(¢) and the 
functions P,(¢) and P.(¢). 

The natural thing to do is to assume N(¢) to be a normal dis- 
tribution. This, however, prevents the evaluation of various integrals 
in finite form. Therefore we shall follow H. D. Landahl (1937) and 
put 


Noo e-7lPl, (9) 


N (6) = 


As Landahl has shown, the results obtained with such a function are 
very similar to those obtained with a normal distribution. 
Equation (9) gives 


Noo et: (10) 


For¢>0 N(¢) =N.(¢) = 
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No oC 
2, 


The function P,(¢) is determined by the distribution function of 
the fluctuation at the connection activated by «, and «,, (Landahl, 
1938). Assuming this distribution function to be 


ere, (11) 


For¢<0 N(¢) =N_(¢) = 


p(§) =F ene, (12) 

we have (Landahl, 1987; Rashevsky, 1948) 
For ¢ > 0 P,=1—4te"?; (13) 
For ¢ <0 yes 7 (14) 


When the situation described by equation (8) takes place then, 
as we have said, we must substitute ¢ + w for ¢. Instead of (18) 
and (14) we then obtain, for y > 0: 


For¢d >—yp Py =1— teow; (15) 
For ¢<—yp P,=4 cow), (16) 


We shall consider ¢, the initial value of the excitation difference, as 
a constant characterizing each individual. Thus P, is a function of 
the initial value ¢ and of the added value w. 

Equation (8) may be written 


Lian Peay eat ome 17 
sts J &@ y($)] d¢—ay. (17) 


But, according to (2) and (3), we now have 
a(d)=Pild,y)N(o);  — -¥(@) =N (6) [1— Pilg, y)]. (18) 


Hence 
«(¢) —y(¢) =N (4) [2P,(¢, y) —1]. (19) 


Introducing now (10), (11), (15) and (16) into the integral in 
(17) we find 


ie [=(9) —v(@)]d6= | N() [2P,($, y) —1] dg 


oO 


c 7 0 
== ‘i ere (ekg) — 1)d¢ as the ert [1 — ekg] de (20) 


: iP e-79 [1 = e* gry) ] ds | E 
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Evaluating the integrals in (20) and substituting them into (17), 
we find, after rearrangements, 


ad. 2 
BLS Hip 1+ e-cy — 


o” — k2 o? — k? 


= c¥} —ay, (21) 


We shall now investigate the expression in braces. There are 
three cases: 


Gasevl. ko ot 


For y = 0 the expression in braces is zero. For y= oo it is equal 
to 1. The second term is negative and decreases to zero more slowly 


FIGURE 1 


than the third, which is positive (Fig. 1, curves II and III). The first 
derivative is equal to 
2 
hk? C pe oy + k oC 
o? — k? o? — k? 
and is positive for y= 0, being equal to ke/(o + k). 
The second derivative is equal to 
k? o? 
o eee ke? 
and is everywhere negative, so that there is no inflection point. 
The sum of the second and third terms in the braces is thus 
shown by curve I in Figure 1. The whole expression in braces is rep- 


resented by curve IV of Figure 1. 


ey (22) 


? 


(eat er?) 
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Case ll. k<o. 


Now the second term is positive and decreases more rapidly than 
the third, which is negative (the curves II and III in Figure 1 merely 
interchange their parts). The first derivative is again positive for 
wy =0. The sum of the two exponentials is thus again a monotoni- 
cally increasing curve, starting at —1 for y = 0 and tending asymp- 
totically to zero (Fig. 1, curve I). The total expression in braces is 
represented by curve IV of Figure 1, and is again the same as in Case 
I. 


Case III. k= ¢. 


Put 
he, (23) 


where A is a small quantity, tending to zero. Preserving only linear 
terms in 4, we have 


f—heS==—2aA. (24) 


Substituting (23) and (24) into the expression in braces of (21) we 
find for that expression, after rearrangements: 


oe genie 
2A 


As A tends to zero, expression (25) tends to 


ad ae (25) 


1— (1+ Jew. (26) 


It is again represented by a curve starting at zero and increasing 
monotonically, tending asymptotically to 1. 


Denoting the expression in braces in (21) by F(w), we see that 
dy/dt = 0 when 


a 
F(y) apne (27) 


If for y = 0 the derivative of AN.F(y) > a, or because of 
(22), if 


a ko 
AN, otk’ 


(28) 


then equation (27) has one positive root *, as is seen from Figure 2. 
The configuration which corresponds to that root is stable, because for 
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yp < y* we have dy*/dt > 0, while for yp>y*,dy/dt<0. 


a 
Ane he 


wr acre WP 
FIGURE 2 


If, however, inequality (27) is not satisfied, then the configura- 
tion X = Y is stable. Though a slight excess of X makes the term 
A(X — Y) in (8) non-vanishing, it remains less than ay. 

To determine the values of X and Y at the stable point y*, we 
must express X in terms of y. The other variable, Y, is given by 


Y=N,—X. (29) 
We have 
X= [- x@as. (30) 


Introducing (10), (11), (18), (14), (15) and (16) into the first 
equation (18), and introducing the latter into (30) we find 


Noo ps4 2 =f ) 
x=—27 14 f or kody + | e9[1—4 eG] dg 
a3 -y 


2 
(31) 
4 jo [1—h eke] a | ; 
0 
Evaluating the integrals, we find 
No o {2 o 1 o i 
x= 5 {2 42 ew Pew 7 | : (32) 


Putting y = y* into (32), we obtain the value of X at that 
point. The problem thus reduces to the solution of the transcenden- 
tal equation (27). The following method of successive approximation 
may be suggested. Equation (27) is of the form 


Ltac%—fpew=yy, (33) 
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with 

k? o ; eee 
iy AN, : 
As a zeroth approximation we take the point of intersection of the 
line y y with the asymptote, that is we put (Fig. 3) 


; a wees (34) 
pe pee oe B pHi] 


FIGURE 3 
1 

yo =. (35) 
i 


Then we introduce this value into F (wy), and take as the first approxi- 
mation the point of intersection of the straight line y y with the hori- 
zontal line F'(y,*) (Line AB, Fig. 3). This gives 


yp=F(y") (36) 
or 
pr — (37) 


Writing F'(y) explicitly, we have 
1 
yy =- (14+ a7 — pe*), (38) 
hs 


For the case in which X is large, this approximation should be fairly 
sufficient. We may, however, proceed further by putting 
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F (y,") 
pt =——, (39) 
? 
and in general 
(wn*) 
ae (40) 


iy 

The larger N,, the smaller y, and therefore, as seen from Fig- 
ure 3, the larger y,*. From equation (32) we see that X increases 
with py. 

We thus obtain the following result: A society in which the 
tendencies toward two mutually exclusive activities are distributed 
amongst the individuals according to a symmetric distribution func- 
tion will, as a result of imitation of behavior, show either one or the 
other of the two behaviors predominantly in certain cases. For con- 
stant neurobiological parameters, A and a, the predominance will be 
greater, the larger the total population N,. For sufficiently large 
values of N, practically all of the society will exhibit the same be- 
havior. The choice, however, between the two mutually exclusive be- 
haviors is determined by pure accident. 

The importance of a possible generalization of this result to n 
different behaviors is obvious. This shall be discussed in a subse- 
quent paper. 

Consider now the case in which the function N(¢) is not sym- 
metric. Let the asymmetry favor the negative values of ¢, in other 
words favor behavior R.. One possible way of obtaining such a dis- 
tribution would be to put, instead of (10) and (11): 


0 01 G2 


N. 
For¢>0 N(¢)=N.(¢) = ey Gas (41) 
No C1 Oo 
For ¢ <0 N(s) =N-(¢) =e (42) 
with 
o1<\ 02. (43) 


Now introducing (41), (42), (15) and (16) into (17) we find 
after elaborate calculations: 
d wy | 2h? C2 


1+ ————c“™ e-ow 
(a+ os) (ox? — kK?) 
(44) 
Steitis satiate ey 
(o, — k) (02 + k) 
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Denoting the expression in braces by F,(y) now, let us investigate 
that function. Let us compare it with the function F'(w) for the case 
in which o = o, in the latter. 

First consider the case in which 


C2 > C1 = k ° (45) 
A comparison of the expressions in braces in equations (21) and 
(44) shows that the second, positive, term in F,(y) is greater than 
the corresponding term in F'(y). We have 
2h? J2 ke k Cg C1 


ee ee eee >0. (46 
(04 at a2) (04? mie k?) oy” PSs k? O71 =F C2 Gc == ke ( ) 


Similarly we find that the third, negative, term of F', (yw) is great- 
er in absolute value than the corresponding terms of F'(w). We have 


0; O02 Cis Oz k Og” OZ 


Ce eee 2 


(04 ca iad k) (2 Sr k) Gas Sear k Len) ait k Cia Pry k2 


But the difference between the absolute values of the positive 
terms is less than the difference between the absolute values of the 
negative terms. Those two differences are given respectively by the 
right-hand sides of equations (46) and (47). Since, because of (45), 


k?<o,k; Oi 1 tei> ost ie, (48) 


a v. (47) 


therefore 


ke O71 k 
< , 
Oz ak C2 G2 Ge k 
and hence the right-hand side of equation (46) is less than the right- 
hand side of (47). 


Thus the negative term in F,(w) exceeds the negative term in 
F'(y) more than the positive term of F,(w) exceeds the positive term 
of F(w). Therefore . 


(49) 


Fi(y) < F(y). (50) 
When 
2 > k = O1,; (51) 
we have 
O71 k k2 
< 


kit+ C02 Cz re: (52) 


Remembering that the second terms now are negative and the 
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third term positive, and that ine iti 
, qualities (46) and (47) are now re- 
versed, we see that the deficiency of the negative term in F, (w) over 
that of F'(y) is greater than the deficiency of the positive term. Hence 
again we have inequality (50). im 
The case k > o2 > o, leads to the same result. In. i 
holds also for y= 0. Hence cares 


POY <0. (53) 


Sa ne y = © do both F(yw) and F,(w) tend to the same asymptotic 
value 1. 


FIGURE 4 


The relations between F (wy), F:(w) and ae 


w are shown graphi- 
0 


cally in Figure 4. If o, become sufficiently small, that is the asymmetry 
favoring negative values of ¢ becomes sufficiently large, other para- 
meters being kept constant, then F,(y) will vary as represented by 
the broken line. It will not intersect the line (a/AN.)y, and a stable 
configuration in which X > Y will not exist. It is, however, readily 
seen that a stable configuration with Y > X will exist in that case. 
Let N(¢) be originally symmetric, (o, = o2 = a), and let the 
stable configuration with a strong preponderance of X , (R:) be acci- 
dentally established. Then let, due to a biological variation of the so- 
cial group, N(¢) become asymmetric in favor of R.(o, < o.). In 
that case as o. increases it will eventually reach the threshold, at 
which the broken curve in Figure 4 will just become tangent to the 
straight line. At this moment dy/dt will become negative and the 
configuration will “jump” over into one for which R, prevails and 


Ve eae 
Another way to introduce an asymmetric distribution function 


N(¢) is to put 
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Noo 


N(¢) = el bo-9| (54) 
The function N(¢) is now symmetric with respect to ¢ =¢., but not 
with respect to zero. Whereas in the case of the function (9) the 
average value ¢ of ¢ is zero, in the case of (54) it is a function of ¢o. 
By the same procedure as above we may derive X and Y as a func- 
tion of ¢ or of ¢o. 

From (18) we see that x(¢) is also a function of the parameter 
wp, and should be written x(¢,y). Since ¢ is independent of ¢, we 
have 


da(p,y) _ eu(p,y) dp 


, 55 
dt op dt Co 
Substituting for dy/dt its expression from equation (8) we find 
dx ( ’ ) ox ( , ) Ox ; 
GND PN a Ey fg 
dt Ow op 


Integrating with respect to ¢ between —o and +o we find, because 
of (380): 


dX kee Me es 
PIE Te a EE genes GENS Hi ae (57) 
dt Es oy : oy 


oo 


We shall now compute the quantity 


+00 ace ( 
Z= 4, AEE oe (58) 
Logt teh 
Introducing (15) and (16) into (18), we find 

x(¢,p) =N (4) [1—$ e*O] for >—y; (59) 

x(d,p) =4N(d) OM for¢<—y. (60) 
Differentiating we have: 
ox(d,y) ik 

i eee =5 N(e) ekou forg >—y; (61) 
ax(d,y) k 

ao Ew) eo forg<—y. (62) 


Comparing (61) and (62) with (59) and (60) we see that 
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ox (d, p 
For $>—y SY SIN) — het, 3 (63) 
Ox (¢, yp) 
For ¢<—y ae a) (64) 
Hence 


—W +00 +00 
Z=k [ a(p,ydetk [ N(s)dgo—k [ x(s,p)d@. (65) 
2 -y - 


Introducing now for N(¢) the expressions (10) and (11), using 
(18), (15) and (16) and evaluating the integrals we find: 


_Noke { oad 


Z, 
4 |e—k o? — i 


eo} : (66) 


Eliminating y from equations (32) and (66), we would obtain 
Z as a function of X. Unfortunately it is impossible to eliminate yw 
because both expressions are transcendental. We may however per- 
form this for the two limiting cases, namely, when o << k&, and when 
ae ae 

When o << k& we have from (66) and (82): 


ok? 
a Ss n eo: (67) 
2 (k? — o?) 
Noo{ 2 o 1 | NK 
— - —- ) er! |=N, + ———— ce. (68 
ap Seer ak 0 3 (@—B) ey 
Hence 
Z=a0(No—X). (69) 
For oc >> k we find in a similar way 
Z=k(Ns= xX): (70) 
Denoting by 1/x an appropriate average of 1/k and 1/o, for in- 
stance 
==" (7 +=) ; (71) 
kin 2 Ke, <0: 


we have approximately, in general, 
Z=x«(No.—X). (72) 


When o << k, we have from (32) 
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1 0 
=-lo —-, (73) 
ee (NGG) 
When o >> k, we have 
i No 
= - log —-_————_—.. (74) 
Ye eee (Meee) 
Hence approximately we have, in general, 
No 
= - log ——_——_., (75) 
omaha CN EEC) 
Introducing (72) and (75) into (57) we obtain 
dX 0 
—=Ax(N.—X) (X—Y) —a(N.—X) log ———_———.._ (76) 
it «(No — X) ( ) —a(No— X) eo. X) 


Because of the assumption of y > 0 in (15) and (16) this equation as 
well as all preceding ones hold only in the range between X = 4 No 
and X =N,. 

For X <4N.; Y > 4 we have, because of the symmetry of N (¢) 


dY No 
—=Ax(No— Y)(Y —X) —a(N, — Y) log ———-,, (77) 
ap ane ) ( ) —a( ) lon ere 
which, because of Y = N,— X , gives 
dX No 
—=Ax«X(X—Y) + aX log —. 78 
dt « X ( )+a cae (78) 


We now may generalize our results by considering not merely 
the frequency of reactions R, and R; , but also their intensities. 

For the two categories case (Landahl, 1938) considered here, it is 
seen that for a constant value of ¢(> 0) the intensity of R, is pro- 
portional to ¢, so that 


R,=a¢. (79) 
A fluctuation — & changes ¢ to ¢ — &, and therefore we have 
Hr a(h =e), (80) 
To have any reaction R, at all; we must have 
aa 60) or E<¢. (81) 


Hence if p(&) is the fluctuation function [equation (12]), we have for 
the average intensity R, of R,: 
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t°) 
a f (@—£)p(é)dé 
er ce ae a (82) 
| pede 


Substituting for p(é) its expression as given by (12), remem- 
bering that 


k 
foré<0 p(é) =e (83) 


k 
fon e.>.0 p(é) ape (84) 


and evaluating the integrals, we find 


2a¢ eee 
F > 22 SSS 
or¢>0 R= — (35) 
sae a 


The average intensity of R, increases with ¢ for positive values 
of ¢, and for very large ¢ we have 


Pods (87) 


For ¢ < 0, when R, is always “accidental” its average intensity is 
constant and equal to that for ¢ = 0, as obtained from (85). 
Equation (82) may be written: 


»p i) 
ad | p(é) dé—a [ Ep(é) dé 
ee a (o— 2). (88) 


fewas 


where ~ denotes the average fluctuation in the interval — o to ¢. 


Computing & , we find 


(89) 
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= i 
For ¢ < 0 el ger (90) 


Substituting (89) or (90) into (88) we again find (85) or (86). 

Now we may consider that the stimulus which increases y, ac- 
cording to equation (8), is not X — Y, but the difference of the total 
intensities R, and R, , that is 


[RG +v)26+v)-—R@t+ y+ wldy. OD 
As is readily seen, this expression is a function of yw only. 

Denote by n the average number of acts of either R, or R, done 
by an individual per unit time. If the average age of an individual 
in the society is 7, then the average number of acts performed by the 
individual during his lifetime is nT. If the probability P, for some 
individuals is less than 1/n7, those individuals will on the average 
never perform a reaction R, in their lifetimes. Similarly, if P, <1/nT, 
such individuals never do, for all practical purposes, perform R, . Let 
A, be the root of the equation 


1 
EQ 8 Be) (92) 
nT 
when solved for ¢ while A, is the root of 
1 
Baht Vs) gana (93) 
nT 


when solved for ¢. Then we may say that all individuals with ¢ < A, 
will never perform R,, while all individuals with ¢ > A, will never 
perform R, . 

We may now break the range of integration in (91) into three 
intervals 


send ioe ee Fi A, , As; As, 80 < (94) 


Denote by a and a appropriate average values of R,(¢ + y) in 
the intervals A,, + o and A,, A, correspondingly. Similarly denote 


by ¢ and ¢ appropriate averages of R, (¢ + w) in the intervals — o , 
A, and A,, A, correspondingly. Putting 


(ce) Ay 
X= | “(o+y) dd; Y= f vty) as; 


ie “ (95) 
x= iE a(p+y) do; y= J y(o+y) de, 
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and using the mean value theorem, we see that the integral (91) be- 
comes 


TOE a ax aaa Cals aa CY . (96) 


Instead of equations (76) and (78) we now obtain: 
For X >4N, 


dx 
yee oe) (G.X) + aX — G)Y,.— éY) 
(97) 
Ne log ————___—_ > 
o—X) eo. X)' 
For X < iN, 
dx 2 No 
GE 8 & (Geko F aX — mY. — c¥) OSS (98) 


If a is so small that the log term may be neglected, and if, as a 
very crude approximation, we substitute in both (97) and (98) an 
appropriate constant time average for (N — X) and X, then both 
(97) and (98) reduce formally to equation (2) of chapter iii of our 
book (Rashevsky, 1948a). As has been remarked in a previous paper 
(Rashevsky, 1949), the coefficients a), ¢, a and c are still functions 
of time, and appropriate averages of them should be used. 


LITERATURE 


Landahl, H. D. 1938. “Contributions to the Mathematical Biophysics of Error 
Elimination.” Psychometrika, 3, 169. 

Rashevsky, N. 1948a. Mathematical Theory of Human Relations. Bloomington 
(Indiana): The Principia Press. 

Rashevsky, N. 1948b. Mathematical Biophysics. Revised Edition. Chicago: Uni- 
versity of Chicago Press. 

Rashevsky, N. 1949. “Mathematical Biology of Social Behavior; II.” Bull. Math. 
Biophysics, 11, 157-163. 


‘ae cor ete 


gree 
Sustingn atte a aaeemaile 
7 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 11, 1949 


OUTLINE OF A PROBABILISTIC APPROACH TO 
ANIMAL SOCIOLOGY: II 


ANATOL RAPOPORT 


COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


Under certain assumptions concerning the probabilities of “muta- 
tions,” ie. changes of structure of bird societies, it is shown that the 
probability distribution for all possible structures of a society of N in- 
dividuals approaches a limit independent of the initial probability dis- 
tribution. A formula for the limiting distribution is derived. 


In a previous paper (Rapoport, 1949), hereafter referred to as 
I, we examined the possible “structures” of bird societies as deter- 
mined by the distribution of peck rights among the individuals of a 
flock. It was shown that on the basis of certain assumptions concern- 
ing the probabilities of the outcomes of encounters between individ- 
uals, it was possible to compute the probability of establishing each 
of the various structures. It was assumed in J that as a result of the 
first encounter between a pair of individuals, the peck right relation 
was established for this pair once for all. 

In the present paper more general assumptions will be made. A 
peck right relation between two individuals will be considered undis- 
puted during a finite time only. Eventually another encounter may 
occur, as a result of which the peck right relation may be reversed. 

Consider a society of N individuals. Associated with N there are 
n possible structures as described in I. Let us suppose that from time 
to time encounters take place between pairs of individuals. The re- 
sult of each encounter may be either the preservation of the old peck 
right relation between the two individuals or its reversal. If the peck 
right relation is preserved, then certainly the structure of the society 
is also preserved. On the other hand, if the peck right relation is re- 
versed, the structure may be changed, or it may not. For example, 
the structure S,: (2, 1, 0) of the three-individual society is diagra- 
matically represented as follows (Cf. I): 
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2 


Of ee 


FIGURE | 


A reversal of peck right between individuals 2 and 0 changes the 
structure to S.: (1, 1, 1), thus: 


<——_ 


FIGURE 2 


; But the reversal of peck right between 2 and 1 results again in 
the structure (2, 1, 0) with the individuals simply relabeled: 


© 


FIGURE 3 


We shall refer to changes in structure as “structure mutations” 
or simply as “mutations” and will denote the mutation S; > S; by 
Si;. The probability of the occurrence of Si; will be denoted by ai;. 

In general, not every mutation S;; can be accomplished by a sin- 
gle reversal of peck right. The corresponding a;; will then be zero. 
Otherwise a;;(i # 7) will depend on the probability of an encounter 
between a pair of individuals which may result in S;; and on the 
probability of the victory going to the “submissive” individual. How- 
ever, in estimating the ultimate fate of the society (the limiting prob- 
ability distribution of its possible structures), one may use the para- 
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meters ai; directly and compute the limiting distribution in terms of 
these. The computation of the a; ; 1n terms of the probabilities of en- 
counters and reversals is a separate problem. We shall first state the 
problem of determining the limiting distribution in terms of the Qi; 


The Matrix (aj;). 


Let Si(t) be the probability of occurrence of the structure Si 
at the time ¢. Take as a unit of time the average interval between 
encounters. Then 


S;(t) = 4118; (¢—1) + oiS,(t—1) + ---- GniSn (t — 1) 
S2(t) = @y2S, (6 —1) + aS, (t —1) + +++ AnoSn(t—1) 


(1) 


UE ay BAS 2 oat [NaS Cael ae ae Gecsn(t——1)s 


Note that the S;; are the “identity mutations,” thus, ai; is the 
probability of the preservation of the structure S;. We denote the 
matrix of the probabilities by (ai;). The attention of the reader is 
called to the fact that the subscript notation is the reverse of the con- 
ventional one, a;; being the element of the 7th column and the jth 
row. We have chosen this departure from convention in order to keep 
the suggestion that a;; stands for the probability of S; mutating to 
re 

If S(t) is the vector {S,(t), S.(¢), ---- S,(t)}, equation (1) may 
be written in vector-matrix notation thus: 


S(t) = (ai;)S(é— 1). (2) 
Hence by iteration, 
S(t) = (aij) S(O), (3) 
and 
S(o) = Lim (ai;)'S(0): (4) 


If the limit in equation (4) exists, then the vector S(co) will 
represent the distribution of the structures S; after a long time, that 
is, if a great many societies are allowed to exist for a sufficiently long 
time, the components of Sco), Si(co) represent the expected frac- 
tions of the total number of societies to have the structure S;. 
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Let us examine the properties of the matrix (ai;). As was shown 
in I, each structure is obtainable from the “simple chain” by a finite 
number of peck right reversals. Therefore the simple chain is so ob- 
tainable from any structure, and hence any structure is obtainable 
from any other by a finite number of reversals. Thus, any change of 
structure S; > S;, even though it may not be possible by a single 
mutation S;;, can take place by a finite number of steps, 


Si Spm ses Sj. 


This is stated in Property 4 below. We shall also assume that any 
structure may be preserved in an encounter. This is equivalent to the 
assumption that the diagonal elements of (a;;) are not zero. Finally, 
since the elements of each column exhaust all possible mutations of 
a structure, all the column sums must be unity. We summarize these 
properties as follows. 


Property 1; QS¢; 21. 
Property 2: Saj;=1. 
j 


Property 3: Waj;;#0. 


Property 4: If ai; = 0, there exist elements @ix, Qx1, Qim,-°* 
@»; , Such that the product @ix@xi +--+ Ami ~O. 


Our principle result is embodied in the following 


Theorem. There exists a matrix (ai;) = Lim (a;;)t, such that 
t=co 
the columns of (ai;) are all identical. The vector represented by any 
of these columns is the limiting distribution vector S(«), and it is 
independent of the initial distribution S(0). 


Proof of the Limiting Distribution Theorem. 


We shall establish a number of lemmas, which will enable us to 
prove the limiting distribution theorem. 


Lemma 1. The elements ai;“ of the matrix (ai;)* satisfy Prop- 
erty 1 and Property 2 of (ai;). 


We shall prove an induction on t. Assume that the elements 
ai; satisfy both properties. 


But Sat int is: : Piet ji the ae 
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.(@) — - = 
ai =o ye at Qe?) Aig + -+-- Oni Ain 


aE) = * 
Bia = Aye Gig + Agg(t-D Dig te9*+ Ane Qin 


: : (5) 


Since 0 = ai; and 0 S aj;*, obviously 0 S ai,“ 


. To obtain 
> i; = 1, we add the equations (5) and get 
i 


> ij = iy > yj" + Ais Dd aaj 
; 


is) J 
(6) 
“aes Qin d cs => ij = 1. 
j j 
This establishes the induction and proves the lemma. 
Lemma 2. For all t , Max. {am , «++» aon‘ } S Max, {an , --+- 
ann} and Min. {ay , ---- am} 2 Min. {ay , ---+ amt}, 


Bea (1. Oy -nosn)s 


Proof. To fix ideas, consider the elements of the first row of the 
matrix (ai:;)*. We have 


Oar = Age ay + dey ys + e+ ++ Oar Ain 


Oey ©) = Ag gy 1 doy Aes + oe 4+ On) Aan 


: (7) 


Bg Og) Any ee oe a ee Ani) Ann » 
Let Gq) = Max. aj‘). Then, by (7) 


Oy) S Aa > ay = A , (8) 
g 


since Sa,;=1. Similarly a, = at, etc., for allan (@=1,2 
; i i @m< (1) = M (#-1) 
---» ), and, in particular, Max. ai = Qa == Max.a,,‘*. 
The proof for the minimum element is exactly analogous. 


Lemma 3. The maximum (minimum) element of each row of the 
matrix (a,;)' tends to a limit as t increases without bound. 
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This follows from the fact that the sequences 
Max. (ay; 5 doi «+++ Gy} and Min. {ai , doi ++++ dns} 
are both monotone and bounded by zero below and by unity above. 


Lemma 4. There exists an integer r such that a,;“ # 0 for 
t> reGoj— 1, 23. - n). 


The lemma implies that sufficiently high powers of (ai;) have 
no zero elements. 


Proof. To begin with, if a;; #0, then ai; #0. We show this 
by induction, assuming a;;“"") #0. Then 


ais = A; Oia + Aq; (F-?) Qis if... ait? ii 
(9) 
5 ey hed) Qin « 
But ai; # 0 by the hypothesis of the induction, and ai; # 0 by 
Property 3. Furthermore, all the other terms on the right side of 
(9) are non-negative by Lemma 1. Therefore ai; # 0. 

Suppose now that a;; — 0. Then by Property 4, there exist ele- 
ments Giz, At,» Gm, **** Gmj, all positive. In conventional notation, 
these elements would be denoted by Qjm, Qmz, Qix, -**: Gri. Now 
aj? =D Andy A 0, SINCE AjmOnmi ~ O. Similarly aj, =D aj, an F 

h h 


0, since aj, a,x # 0. Proceeding in this way, we obtain a;;* # 0 
for some integer 7. But this is a;;“ in our notation. It follows by 
the proof of the first part of the lemma that all higher orders a;;“ 
are not zero. This proves the lemma. 


Lemma 5. Let a; = Lim. Max. {a,;", ao; ---- a,j}. Then 


t=co 
for every « > 0, there exists an integer t,, such that fort > t, 
ja; — ai?| < e (i =1,2----n), that is, all the elements of the jth 
row of (ai;)* tend to a; as a limit. 


Proof. Again let us fix our attention on the first row. Let s be 
such that for t > s Max. {ai} — ao, < e,, where «, =e - Min. 
aij/(1 + Min. ai;™) as shown in Lemma 8. Let 7 be such that 
for t 2 7, ai; > 0, as shown in Lemma 4. Choose t = r + s. Then 


(ais)* = (aij)* (Qij)". Let ae = Max. {d4,, doi «+++ Gai} and 
Qin = Max. {du, de® +--+ Qu}. Then 
Cie = Air Ama" Fier? Ame? ++ Qa Ann”, (10) 


since + @m; = 1 by Lemma 1. On the other hand 
j 


Qa = Ary Ama" + Aer) Ame’ + +++ + Any Ginn’. (11) 
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By Lemma 2, «; S dni S ay,, so that by our hypothesis on s , 
is) — Am < e,. Subtracting equation (11) from (10), we get 


Omi) (A? — day) + One (dx) — 23 ‘*)) 

+ .... Qn’? (dy = Oni) < oh (12) 
We have proved the first part of our theorem, namely that (aji;)¢ 

tends to a limiting matrix with equal columns. For sufficiently large 

t, therefore, (ai;)* may be written as (a:i;) + (e€;) where 


Oy Oy ++ Oy 
Oz Og +++: 

iy = (13) 
An eeeeee Qn 


and |Max «;;| is arbitrarily small. Now if S is an arbitrary vector 
with non-negative components, whose sum is unity, the product 
(ai;)S is independent of S, being the vector {a} = {a,, O-:-- Gn}. 
On the other hand («:;)S gives the vector {> e811, 5 eieSe, ---- 


7 v 
> éinSn}', each of whose components is less than |Max. «;;| in abso- 


lute value and hence arbitrarily small. But (ai;)*S = (ai;)S + (e;)S 
= {a} + infinitesimal. Therefore Lim (ai;)‘S = {a} and is independ- 
t=00 


ent of S. Our theorem is thus proved. 


The Calculation of the Limiting Distribution. 


We have shown that for every mutation probability matrix hav- 
ing properties 1 - 4, there exists uniquely a limiting distribution 
vector. This vector is invariant under the linear transformation 
(a;;) and therefore satisfies the equations 


(ai))S=S; (14) 
DS:=1. (15) 


The solution of (14) reduces to the solution of n — 1 simultane- 
ous linear equations. The existence and the uniqueness of the solu- 
tion is guaranteed by the limiting distribution theorem. We have, in 


view of (15) 
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Si = Q1uS1 + GaSe + ++++ Ani (1—> Si) 


i=1 


SY] US abST aS (16) 


4-1 


n—1 
Sa = Qi (m1) 91 a G2 (n-1)S2 + ++ On (na) (1—3S;) ? 


q=1 
and, after proper transpositions, 
S,[1 — (day — Gn) | + So (Qe1 — Gq) +---- 
Si (ai Qin-1)1) = Ong 


S2(@ne— Giz) + So[1 — Ges — Ane) ] +--+ 
Sna (Anz — @in-ay2 ) = One (17) 


Sort) (Gn (n-1) — G1 (n-1) ) tosses 
S(n-1) [1 me (Gn(n-1) Tz Qin-1)(n-1)) J = Qn(n-1) « 


Denote by (bi;) the (n — 1)-rowed matrix obtained by deleting 
the nth row and mth column of the matrix, [I — (ai; — Gnj)], where 
I is the identity matrix and denote by (b;;)‘ the matrix obtained by 
replacing the ith column of (bi;) by the vector {@n1, Gyo, «+: Qn (n-1) }- 
Then by Cramer’s rule, the components of the limiting distribution 
vector are given by 


| (bis) ®| 
So 


, (4=1,2..--n). (18) 
| (bi;) | 


Example: Society of Three Individuals. 


In case N = 8, n = 2, and the system (16) reduces to a single 
equation, 


Si = ayS, ahs G21 (1 — S,) (19) 
whose solution is 


Si = Qei/(1— ay + ex), (20) 
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then obviously 
S,=1—S,= (l—a@n)/(l1—au + Qe1). 


We can now make assumptions concerning the probabilities of 
encounters and victories, and on the basis of these assumptions com- 
pute a,, and ad ,. If an encounter between any pair of individuals is 
equally likely and the probability of victory does not depend on the 
peck right relation existing before the encounter, that is, it is 4 for 
each individual, the calculation of the aj; is quite simple. Since the 
mutation S., can occur as a result of any encounter (Cf. Fig. 2), 
provided the peck right relation is reversed, we have a,, = 4. On 
the other hand, S, is preserved in 5 out of 6 possible results of en- 
counters as can be seen from Figure 1. Hence a, = 5/6, and 


= 1/2 = 
et 6 / Ge 17/2 0. 


Note that this distribution is also the initial distribution for N = 3 
shown in I to result from random initial victories. 

The “mutation” method, however, enables us to introduce biases 
which may depend on inherent properties of individuals and on their 
social rank. Thus the probability of an encounter between individ- 
uals of widely different social rank may be taken to be smaller than 
that between individuals of nearly equal rank. Likewise the prob- 
ability of victory of a dominant individual may be taken to be greater 
than that of the submissive individual, etc. Hysteresis phenomena 
may likewise be introduced, that is the dependence of victories on the 
number of past victories enjoyed by the individual, etc. Some of these 
complications will be discussed in later papers. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 


3/4; S,=1/4. (21) 


1 
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Some neural mechanisms are described which interpret neurobio- 
physically the determination of the behavior of an individual by the 
maximizing of his satisfaction, or pleasure. 


The concept of the satisfaction function has proved to be very 
useful in the theory of social behavior (Rashevsky, 1948a; Rapoport, 
1947a,b,c; Rapoport and Shimbel, 1947). The postulate that an in- 
dividual behaves in such a way as to maximize his satisfaction func- 
tion, or the sum of the satisfaction functions of all individuals, has 
been introduced in a purely formal manner without any neurobio- 
physical interpretation. An attempt at such an interpretation has 
been made by H. D. Landahl (1948). However, Landahl’s theory 
has an unpleasant feature, namely, that the behavior of an individual 
is determined not exactly by the maximum of certain neurobiophysi- 
cal functions, but only approximately. In the present paper we shall 
suggest two mechanisms which avoid this difficulty. 

The first mechanism assumes that hedonistic behavior, that is 
behavior characterized by the maximalization of one’s own satisfac- 
tion function, is not innate, but learned. 

Consider the case studied by H. D. Landahl [1941; Rashevsky, 
1948b (chap. xli)] in which an individual learns to choose one “cor- 
rect” stimulus out of a total of N stimuli. In Landahl’s theory the 
“correct” stimulus is characterized by the circumstance that the re- 
action to it produces an additional central excitation which reinforces 
that reaction. A “wrong” stimulus is characterized by the fact that 
a reaction to it results in a central inhibition which weakens the re- 
action. 

Let us now modify that situation somewhat by considering that 
every reaction to any one of the N stimuli results in a reinforcing 
excitation. The degree of reinforcement is, however, assumed to be 
different for the different stimuli. The stimulus which is reinforced 
the most strongly will prevail, gradually inhibiting the other stimuli. 
Thus the organism learns to choose the stimulus which results in the 
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greatest reinforcement. Considering the degree of reinforcement as 
the measure of subjectively felt satisfaction, or pleasure, we find that 
the organism learns to choose the stimulus which results in the great- 
est pleasure. 

The rate of this learning can be obtained directly from Landahl’s 
equation by making the constant 6 [Rashevsky, 1948, p. 473, equation 
(11) ] negative. 

Now let us consider that the same reaction of different intensity 
produces different environmental changes in such a way that the 
change produced by a certain optimal intensity results in a maximum 
pleasure. Since stimuli of different intensities result in differently 
localized central excitations (Rashevsky, 1948b, chap. xxxiii), the re- 
actions to different intensities are produced via different central re- 
gions. A choice between different intensities of the same reaction is 
thus reduced to the choice between spatially different central stimuli. 
The case is reduced to the preceding one. Again the organism learns 
to use that intensity of a reaction which gives maximum satisfaction. 

We may finally consider the case in which a constant stimulus 
produces a constant reaction, resulting in an increase of something 
which itself acts as a rewarding stimulus. This something may have 
a maximum value as a function of the duration of the reaction. 

Let the organism cut off the reaction at different times ¢,, t. 

-- t, , among which is included the optimal time t*. To each of these 
times there corresponds an intensity, or quantity, S,, S. ---- S, of 
the “something” that results in some satisfaction. Thus the organ- 
ism develops a series of conditioned reflexes, in each of which a given 
intensity S; results in the cutting off of the reaction. If all those 
conditioned pathways are cross-inhibited, then the strongest will pre- 
vail. But that strongest conditioned pathway will be the one corre- 
sponding to S*, the strongest conditioned stimulus. 


FIGURE 1 


Now let us consider a different mechanism which does not involve 
learning directly. Let the stimulus S excite the center P, which we 
shall call “pleasure center” (Fig. 1). Let the intensity of excitation 
of Phave a maximum for a certain intensity of S. 
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Let the pleasure center P excite a pathway, I, which leads to 
the connection s, and which is characterized by the following rela- 
tion between its parameters (Rashevsky, 1948b, p. 444). 


Be tA eb <a: eaee (1) 


At the connection s there is a self-circuited neuroelement which pro- 
duces a cumulative increase of excitation, resulting in the reaction 
FR. The latter in its turn again produces the stimulus S. 

Under those conditions any sufficiently sudden increase in the 
excitation of P will result in a temporary excitation of the pathway, 
and therefore in a reinforcement at s (Rashevsky, 1948b, p. 444). 
On the other hand, any sufficiently sudden decrease of the excitation 
of P will result in a temporary inhibition of the pathway, and there- 
fore in a weakening of the cumulative excitation at s. 

Let the organism produce a reaction R resulting in a certain S. 
This will produce a certain P, and as a result R will be reinforced. 
Thus a stronger S will be produced the second time, which will again 
be reinforced as long as an increase in S results in an increase of 
excitation in P. But when the excitation of P reaches a maximum, a 
further increase in R, or S, will result in a decrease of the excita- 
tion at P, and therefore in a decrease of reinforcement at s. Thus 
R and S will be reinforcing each other only until the excitation at P 
has reached a maximum. An accidental increase of R beyond the 
value Rnax Will result in a decrease of excitation at P, and therefore 
in a weakening of R. 

The mechanism requires sufficiently sudden changes in the in- 
tensity of R, or S. It is, however, likely that it may be generalized 
to continuous changes. 

According to the picture suggested here, an organism acts altru- 
istically, that is tries to maximize the satisfaction of others, if the 
satisfaction of others acts as a pleasant stimulus on that individual. 
From this point of view altruism is just another form of egoism. 
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Various types of analyses of the existing data on spontaneous and 
X-ray induced mutations of the tobacco mosaic virus to the aucuba va- 
rlant are discussed on the basis of the known phytopathological, bio- 
chemical and biophysical facts. It is concluded that if virus protein is the 
X-ray absorbing medium through which the mutation is induced 3% of 
the total virus protein would be involved in this primary process. If 
nucleic acid were the medium through which the ionizing energy must be 
absorbed in order to induce the mutation, 25% of the amount of this acid 
existing in the virus would be needed. The analysis shows also the com- 
plete insignificance of the naturally occurring ionizing radiations for the 
process of spontaneous mutation of the virus. 


Since the rediscovery of Mendel’s idea early in the current cen- 
tury the geneticists accepted as their fundamental concepts the cor- 
puscular structure of the hereditary material in the cell and the ran- 
dom and discrete character of the primary processes which confer to 
the nucleus a hereditary change. During the same time physicists 
discovered the corpuscular structure of radiation energy and the ran- 
dom character of its interaction with the irradiated atoms. The dis- 
crete and haphazard character of directly observable events induced 
by the absorption of radiation energy has been revealed in a convinc- 
ing manner by the cloud chamber, the electronic counter and the 
photographic emulsion. With less direct evidence but perhaps with 
equal firmness the concepts of corpuscular structure, of a random 
type of action and of a discrete number of essential primary events 
became basic ideas in the study of hereditary changes. It is therefore 
perhaps not chance that the applications of radiations in genetics 
have been so rich in results; since the concepts of randomness and 
discreteness are fundamental to both radiation physics and genetics. 

Whatever the form of the energy inducing a mutation is and 
whatever may be the mechanism of the mutation the experimental 


facts suggest that: 
(i) in order to induce a mutation a certain amount of energy 
must be absorbed by a certain volume (v) in the cell; 
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(ii) the absorption of this amount of energy by that volume 
consists of a number of particular random events (E). 


It is generally believed that the number of such events which 
are sufficient to induce a mutation is small, mostly one for mutations 
occurring in a single gene, and one or two for changes involving a 
whole chromosome (cf. Timoféeff-Ressovsky, 1937, pp. 122-123 ; Fano 
and Demerec, 1944; Lea, 1947a, 1947b). However examples are avail- 
able in which a larger number of primary events is necessary for the 
induction of a detectable mutation. Recessives in polyploids belong 
fundamentally to this case; in the field of unicellular organisms cases 
are known in which effects essentially equivalent to lethals appear as 
a cumulative action of several events (see e.g. Lacassagne, 1930; 
Glocker, Langendorff and Reuss, 1933; Schreiber, 1934). Consequent- 
ly we will not put a priori any restriction on the number of events 
(E), since this should be obtainable from an analysis of the experi- 
mental data. 

Events (E) physically identical with those inducing a mutation 
occur in the whole organism and throughout the entire sample under 
experimentation but only few of them occur in the volumes (v). We 
will indicate with p the probability that an event (EF) occurs in the 
volume (v) of a given biological unit (a sperm, a cell, a whole or- 
ganism) if 7é occurs in the experimental material containing that unit 
once per each u® in average. The probability » for the induction of 
a hereditary change depends on the geometric structure and on the 
chemical composition of the volume (v) and of the whole experimen- 
tal object. 

We will indicate with N, the number of events (E) produced 
per «? of the experimental material by the amount of energy D. We 
will assume that each of these events, if it occurs in the volume (v), 
has an immediate effect, although its detection may be delayed. To 
take into account the spontaneous mutation we have to keep in mind 
that when the administration of the experimental amount of energy 
D starts, some of those effects may be already existing in some units 
of our sample. We will indicate with Ny, the number of ioniza- 
tions that would have to be produced per u? of the sample, in order 
to produce exactly the effects existing in it, if the sample were com- 
pletely free from such effects. With respect to the type of energy 
actually used in the experiment and with reference to the phenotypic 
character under consideration, Ns is some sort of energy equivalent 
of the genetic status of the sample. We will assume that the sample 
is a statistically correct representative of the population out of which 
it is taken and that this population is in equilibrium with respect to 
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the mutation considered, so that there is always a constant prob- 
ability for a spontaneous mutation. We will indicate with P(Ns + Nz) 
the probability of inducing a mutation with an amount of adminis- 
tered energy equal to D, consequently the probability of a spontane- 
ous mutation will be indicated by P(Ns). The term Ns + Np ex- 
presses here the cooperative action of the experimentally adminis- 
tered and of the spontaneously absorbed energies. In analyzing the 
induction of mutations by X-rays, we assume that the events (E) 
consist of ionizations. Approximately 32.5 eV are required to pro- 
duce a pair of ions in biological material, so that in a round figure 
two ionizations per w® are produced by each roentgen of absorbed 
energy; consequently if D is measured in roentgen, we may take as 
a very crude approximation: Np ~ 2D. 

Two situations at least are possible a priori: (i) the appear- 
ance of a mutation may require the occurrence of an exact number 
n of events (E) in (v); (ii) the genotype or at least its phenotypic 
expression may be stable to such an extent that events (E) occurring 
in (v) beyond a certain number x» do not produce any detectable 
change within experimental limits. In the first case the probability 
of induction of a mutation is: 


in the second case it is: 


Po— > 0i— 1 — >, Di, (1) 


i=n 74=0 


where N is the number of events (E), ie. of ionizations, per we of 
the sample. If N is very large there are two types of asymptotic re- 
lations (cf. e.g. Castelnuovo, 1925, chap. v; Uspensky, 1937, pp. 129- 
136; Cramér, 1946, pp. 198-206): 


Dn = (2aN pq)? exp [— (Np — n)?/(2Npq)]; (2) 


P,= f arnt eels s/n (3) 


with u= (Np — n)/(2Npq)*”, (theorems of De Moivre, Laplace and 
Tchebychef) and 
p, =e” w"/n! with w=Np, (4) 
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n-1 
P,=1—e’D> wi/t! 
4=0 


WwW ee 
=| e* 1 dx/(n—1)!=—y(w,n) 


(Poisson distributions). This last integral is the incomplete gamma 
function for which the symbol y is here used. The validity of for- 
mulas (2) and (3) requires in general that p and q be large, where- 
as in formulas (4) and (5), it is sufficient that p or q be small. In 
all radiogenetical applications of probability it has been found that 
is extremely small; for this reason we will use the formulas (4)— 
(5). 

If more than one event (E) is necessary for the induction of a 
detectable mutation the question arises whether the effect of each 
such event is hereditary or not. We will assume that it is, in the same 
way as a mutation in a single recessive gene is hereditary. However 
this question enters into our analysis only insofar as a cellular divi- 
sion may occur between two events (E) in (v). This may happen 
very infrequently in mutations induced by radiations, because of the 
very short duration of experimentation. The question may be im- 
portant though in interpreting spontaneous mutations. 

We proceed now to a more concrete interpretation of the prob- 
ability ». Let us consider for instance the expression (5). If we 
take the number of ionizations produced per yu? as the measure of 
the energy absorbed, the mean energy absorbed per mutation is n/p. 
This is a known result from the theory of the cumulative Poisson 
distribution [cf. e.g. Opatowski, 1946, formula (2)]. If a mutation 
requires » ionizations in the volume v, it requires ”/v ionizations 
per u®; consequently n/p = n/v and » =v. Therefore if the energy 
producing a mutation is measured as the number of ion pairs that 
it produces per x? of the irradiated material, the probability p is the 
volume v in u? in which this energy must be absorbed. to produce the 
mutation. This conclusion implies the assumption that each event 
(E) occurring in (v) produces its effect. If a certain probability of 
effectiveness is attached to the events (E) the interpretation of p 
must be modified accordingly (cf. e.g. Timoféeff-Ressovsky and Zim- 
mer, 1947, pp. 97-99). 

We analyze here the spontaneous and X-ray induced mutations 
of the tobacco mosaic virus from ordinary type to aucuba type (Gow- 
en, 1941). Because the reproductive mechanism of a virus consists 
of a simple division, it is reasonable to assume, at first at least, that 
any hereditary change in the virus is transmitted to both of its off- 
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springs. We assume also that the hereditary changes with which we 
deal do not affect the rate of reproduction of the virus. Then the 
proportions of viruses of a sample containing certain hereditary prop- 
erties are independent of the process of division, consequently this 
will not appear in the mathematical expressions of those proportions. 

The viruses are detected by counting the number of lesions that 
they produce in plants. For viruses used in solution, the number of 
lesions is proportional to the number of viruses only over a small 
range of their concentration. However since in the experiments of 
J. W. Gowen a single type of technique has been consistently used, 
We assume, at first, that his data are not affected by a concentration 
effect (cf. Bawden, 1943; Doerr, 1939; Stanley, 1938). It will be seen 
later that this assumption is not likely to involve a serious error. 

In using the experimental results we have to keep in mind that 
the X-rays have at least two effects on the virus: mutation and in- 
activation, and that reverse mutations occur also. Since no experi- 
mental technique is known which would prevent a simultaneous oc- 
currence of these effects, several types of analysis are possible. 

(i) In view of the fact that the aucuba mosaic viruses, which 
are formed in the course of experimentation from the tobacco mosaic 
viruses, are subject—after their formation—to a much smaller dosage 
of X-rays, in average, than the original tobacco mosaic viruses, we 
may neglect the reverse and other possible mutations as well as the 
inactivations of those newly formed aucuba mosaic viruses, in a first 
approximation at least. However we will take into account the fact 
that a tobacco mosaic virus must survive a dose D in order to be mu- 
tated by it. The data for the inactivation will be taken from the re- 
sults of J. W. Gowen (1940) which give for the probability of a 
tobacco mosaic virus surviving as such a dose of D roentgen of 1.5 A 
wave length the following relation: 


S(Np) = exp (— 415 X 10° D), (6) 


which is valid for the range of D corresponding to the mutation data 
(see Fig. 1). Exponential survival curve is given also by A. Marshak 
and W. N. Takahashi (1942), although differences in experimental 
conditions make their coefficient of D different from that of J. W. 
Gowen. We take P(N; + Np) as the probability of the mutation of 
a tobacco mosaic virus under the assumption that it survives as such 
the administered dose D of X-rays. Then we have to put: 


P(Ng + Np) = A/[TS(Nd)] , (7) 


where T is the number of the tobacco mosaic lesions obtained from an 
unirradiated sample and A is the number of the aucuba mosaic le- 
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sions obtained after the irradiation of the sample; S(Np) is calcu- 
lated here for each D from (6). To correlate the probability of mu- 
tation with the energy which produces it we calculate P(Np) from 
the following equation 


P(NectiNs): =P Neetihs RGN) RNa 2 (8) 


This is a probabilistic expression of the two ways in which a virus 
may mutate; either spontaneously or in consequence of the energy 
D. From (8) we obtain 


P(Np) = [P(Np + Ns) — P(Nsg)1/[1 —P(Ns)] 


(9) 
= P(Np + Ng) —P(Ns), 


since P(Ns) is negligible with respect to unity. In the present case 
it is 0.00079. Applying formula (5) we obtain two equations for the 
unknown v and n: 


y(UNs, 2) =P(Ns) ; 


the right-hand sides being known from observation. 

The lower part of Figure 1 on the right represents the prob- 
ability of mutation P(N>) as a function of the energy in roentgen 
to which the mutation is due. The values represented by dots have 
been calculated by means of the formulae (7)—(9) using the ex- 
perimental data of J. W. Gowen. In order to change the roentgen of 
the figure into Np, i.e. into number of ionizations per u°®, the ab- 
scissae must be multiplied by a suitable factor. If it is assumed 
that the genetically effective ionization occurs in the whole virus 
this factor may be taken as 1.45; this is a value given by D. E. Lea 
(1947a, p. 8) for dried virus at 1.5 A, which are the experimental 
conditions of the analyzed data. The application of the theory to 
these data has been carried out in the following manner: First of all 
from the equation of spontaneous mutations 


y(vNg , n) = 0.00079 (11) 


the dependence of vNs on n has been numerically tabulated. These 
values of vNg have been substituted into equation (10) and from the 
equations thus obtained, a numerical relation between v and n has 
been obtained for each pair of values of P and N>, calculated from 
the experimental data. If for a certain value of n the fit were per- 
fect, the corresponding values of v calculated for all the pairs of 
(P, No) would be the same. However such a perfect case does not 
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FIGURE 1. Probability P(Np) of mutation of the tobacco mosaic virus to 
aucuba against the energy D in roentgen which produces this mutation. 
Theoretical curves and points computed from the experimental data. 
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VALUES OBTAINED FROM THE THEORETICAL FORMULAE 


Theory I (lower part of the figure to the right) 


n v in w3 a/v in % 

1 4.65 X 10-8 45.7 

2 3.80 10-7 14.3 

3 8.97 X 10-7 4.5 

4 1.42 < 10-6 2.0 

5 1.88 X 10-6 PAE 

6 243) 10-8 a3) 

ff 2.94 X 10-6 ee 
100 2.20 X 10-5 12.0 


Theory II (upper part of the figure to the left) 


n vin we a/v in % 
1 1.09 X 10-8 15.7 
2 1.67 X 10-7 21.8 
3 4,24 * 10-7 24.4 


v is calculated on the basis of an average absorption coefficient of the whole virus 
which is practically equal to the absorption coefficient of the protein component of 
the virus. The calculations imply that each ionization in (v) is effective. If it is 
assumed that ions appear in clusters of 3 in average (Jordan, 1988) so that each 
effective event (E) consists not of one but of three ionizations, the above values 
of v must be multiplied by 3. 


present itself in general. Consequently the standard deviation of all 
the calculated values of v from their arithmetic mean @ has been 
computed. The smallest standard deviation is obtained for n = 4 
with o/v = 0.02. Figure 1 and the table show that the fit is very 
good and that it would not be so with another value of n. The theory 
gives in this way a volume v of 1.42 X 10-* uw’, which is practically 
identical with the effective absorption volume of a vital dominant 
gene in the sex chromosome of Drosophila (1.1 < 10-* u® according to 
J. W. Gowen, 1934; experiments with approximately 1.5 A X-rays). 
It is also about of the order of magnitude of the volume for the sex 
linked lethals in Drosophila: 10 wu? according to U. Fano (1942). J. 
W. Gowen (1940) calculated a volume of 7.5 X 10-° u®, whose ioniza- 
tion produces the inactivation of the tobacco mosaic virus. A. Mar- 
shak and W. N. Takahashi (1942) give a size of this volume between 
4.6 X 10° «? and 42 X 10-° uw. The fact that the volume for inacti- 
vation is several times larger than for mutation agrees with the fact 
that it takes much more energy to mutate a sample of viruses than 
to inactivate it. — 

A very large number of papers has been written on the size of 
the tobacco mosaic virus. Methods of ultrafiltration (Bawden and 


x 
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Pirie, 1937; Smith and MacClement, 1938), electron. microscopy 
(Kausche, 1939; Kausche, Pfankuch and Ruska, 1939; Anderson and 
Stanley, 1941; Stanley and Anderson, 1941; Stanley, 1943; Williams 
and Wyckoff, 1945), X-ray diffraction (Bawden, Pirie, Bernal and 
Fankuchen, 1936; Bernal and Fankuchen, 1941), measurements of 
sedimentation, diffusion constants and viscosity (Bechhold and Schles- 
inger, 1933; Wyckoff, Biscoe and Stanley, 1987; Frampton and Neu- 
rath, 1938; Lauffer, 1938a, b; Mehl, 1938; Neurath and Saum, 1938) 
have been all used to find some information on the shape and the 
size of this virus directly or through its molecular weight determi- 
nation (Bernal and Fankuchen, 1937; Bernal, 1938; Stanley, 1938, 
1939, 1941; Stanley and Loring, 1938; Frampton, 1939; Lauffer and 
Stanley, 1939; McFarlane, 1939; Cohen and Stanley, 1942; Mark- 
ham, Smith and Lea, 1942; Lauffer, 1944; Harrow, 1946, pp. 80-84; 
Lea, 1947a, p. 186). The extreme results of most of these determi- 
nations differ from each other by a factor of few units. The differences 
are due, in part at least, to experimental techniques which may cause 
such effects as breakages of the virus by mechanical stress or by 
thermal agitation (cf. Sigurgeirsson and Stanley, 1947). Reasons 
have been also suggested to doubt whether the size of the tobacco 
mosaic virus can be considered a constant (Pirie, 1945, p. 15; 1946), 
since there is some evidence that the virus consists of basic particles 
which arrange themselves in a linear fashion, the number of such 
particles of the virus and consequently its length being dependent on 
the environment (Bawden and Pirie, 1945). However the distribu- 
tion of this length has a pronounced and extremely sharp maximum 
at 280 mu (Sigurgeirsson and Stanley, 1947; Oster and Stanley, 1946; 
Rawlins, Roberts and Utech, 1946). This is also the minimum length 
at which the virus shows its infective activity. It is therefore justi- 
fied to consider it as the length of the virus (Stanley et al). The 
average value accepted today by most workers as the volume of the 
tobacco mosaic virus lies between about 4 X 10° to about 5 X 10° we. 
If the above quoted length of 280 mu together with a cross-sectional 
diameter of 15 mu according to W. M. Stanley et al is accepted, the 
volume comes out to 4.96 X 10-* u*. It is for this reason that we will 
assume in the following calculations the round figure of 5 X 10° yu 
as the volume of the virus. The fact that in fresh unaltered plant 
juice the virus may appear in aggregates of a much larger size (Wil- 
liams and Steere, 1949) has no bearing here. ; 

It has been recognized that a mutation of the virus is due neither 
to a splitting of its original volume nor to a polymerization (cf. e.g. 
Gowen, 1945). The problem arises then to identify within the virus 
the volume (v) responsible for its mutation. J. D. Bernal and I. Fan- 
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kuchen (1941) found through an extensive X-ray investigation that 
the protein of the tobacco mosaic virus consists of many hexago- 
nal unit cells with a base of 87A side length and 68 A height. This 
gives for the volume of each such cell 1.34 X 10° u’, which is very 
close to our theoretical volume of v = 1.42 X 10° uw. According to 
J. D. Bernal ‘and I. Fankuchen the unit cells are arranged along the 
length of the virus with their hexagonal bases adjacent to each other. 
This type of structure is confirmed by recent high magnification elec- 
tron micrographs (Markham, 1947; Wyckoff, 1947). If the ‘“‘muta- 
tion volume” (v) is actually one of the hexagonal cell units the pres- 
ent theory would indicate that these units, although geometrically 
identical, have different biological functions because only a particu- 
lar one of these units would be responsible for the mutation of the 
virus. Such a situation is not surprising in principle, since one can 
imagine physiochemical differences between the end and the inter- 
mediate cells. 

The calculation of the volume within which the primary process 
of mutation takes place gives us the possibility to consider a contro- 
versial question. E. Pfankuch, G. A. Kausche and H. Stubbe sug- 
gested in 1940 that mutation of the tobacco mosaic virus under the 
action of X- and y-rays may be due to changes in its nucleic acid con- 
tent. Their suggestion was based on solubility determinations and 
on other physicochemical measurements, but mainly on electrophore- 
sis which showed that at least 90% of the protein of the tobacco mo- 
saic virus and of its variant strain were identical. The tobacco mo- 
saic virus is known to be a nucleoprotein and its nucleic acid group 
to be particulary efficient in responding to the action of X-rays be- 
cause of its much higher average atomic number due to a relatively 
high content of phosphorus. (The mass absorption coefficient for 
X-rays is proportional to more than the 4th power of the atomic 
number. Nucleic acid of the tobacco mosaic virus contains about 9% 
of 1;P whereas the element of highest atomic number in the protein 
is 169 whose content however is only 0.24%; see Loring, 1939; Stan- 
ley, 1939.) 

On the contrary, W. M. Stanley (1941), C. A. Knight and M. A. 
Lauffer (Knight and Stanley, 1941; Stanley and Knight, 1941; Knight 
and Lauffer, 1942) concluded after their analyses that the differ- 
ences between strains of tobacco mosaic virus are located probably 
in amino acids and that there is little evidence at least for their be- 
ing in the nucleic acid. The question is of interest because it has been 
conjectured that the nucleic acid plays an important role in the in- 
duction of mutations in general. In fact, the dependence of the ab- 
sorption coefficient of the nucleic acid on the wave length of ultra- 
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violet and the dependence of the genetic effect on the same wave 
length show a partial but quite close correlation (Catcheside, 1948, 
p. 358). Even the lethal spectrum of certain unicellular organisms 
has been found to coincide with the absorption spectrum of the 
nucleic acid in the ultraviolet region (Jordan, 1939, p. 479). 

Assuming for a moment that the whole volume (v) which is the 
site of the primary mutational events consists entirely of nucleic acid 
the size of (v) could be calculated in the following manner. The X-ray 
mass absorption coefficient of this acid due to photoelectric effect and 
to transfer of energy to recoil electrons is 12.46 gr> for the wave 
length of 1.5 A used by J. W. Gowen. This value has been calculated 
according to the method indicated by D. E. Lea (1947a, p. 349), i.e. 
from the experimental values of S. J. M. Allen for the total mass ab- 
sorption coefficients (Compton and Allison, 1946, p. 802) the absorp- 
tion coefficient due to scattered radiation has been subtracted. The 
latter has been calculated from a table of D. E. Lea (1947a, p. 348) 
based on Klein-Nishina formula (cf. Mayneord, 1940). The chemical 
composition used for these calculations was H = 4.16%, C = 33%, 
N = 15.1%, P = 9.13%, 0 = 38.61% (by difference). These are av- 
erages calculated from the data of H. S. Loring (1939). On the ba- 
sis of the calculated absorption coefficient of 12.46 gr and of the 
value of 1.6 for the density of the nucleic acid of the tobacco mosaic 
virus (Stanley and Loring, 1938) it is found that 1 roentgen pro- 
duces 3.74 ionizations in 1 «* of the nucleic acid. Consequently if 
that part of the administered X-ray energy which is responsible for 
the mutation were absorbed entirely by the nucleic acid, the volume 
v would be: 


1.42 X 1.45 X 10-°/3.74 = 5.5 X 107 p, 


because in applying formula (10) to a set of experimental data, a 
change of the number of ionizations per u? per roentgen implies a 
change in the same proportion of Np and consequently a change of 
vy in the inverse proportion. 

Since the nucleic acid amounts to about 5.5% in weight of the 
virus (Cohen and Stanley, 1942; also Schramm and Baunendere 
1944) and the density of the latter is 1.3 (Stanley and Loring, 1938; 
Lauffer, 1944) the volume of the nucleic acid existing in the virus 1s, 
if the size of the latter is taken as 5 X 10° y° 


5 X 10° X 0.055 X 1.3/1.6 = 2.23 X 10° p’, 


ich i i 4 ired for the mu- 
which is about 4 times larger than the volume y require 
tation by energy absorption considerations, under the assumption 
that this volume consists entirely of nucleic acid. We may conclude by 
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saying that if the nucleic acid is responsible for the mutation of the 
tobacco mosaic virus, a substantial part of this acid existing in the 
virus, about one fourth, must be active in absorbing the energy needed 
for mutation. 

However if the protein were the only site of the primary muta- 
tional process, the volume involved would amount to 


1.46 X 10°*/(5 X 10% — 2.23 X 10°) = 3% 


of the whole virus protein which would be still in agreement with 
the measurements of E. Pfankuch (1940) who concluded that at most 
10% of the protein may be different in the tobacco mosaic virus and 
one of its variant, although the reliability of Pfankuch’s electropho- 
retic method is put in doubt by C. A. Knight and M. A. Lauffer 
(1942). 

It is seen in this way that the present analysis is not in disagree- 
ment with either of the two alternative theories as to the site of the 
primary mutational process (Stanley et al and Pfankuch et al). It 
may also be that both nucleic acid and the protein are involved in 
the mutation of the virus. W. M. Stanley and C. A. Knight (1941) 
suggested that the mutation of the virus may be a cumulative effect 
of events occurring during its multiplication and it is known from 
higher organisms that both nucleic acid and the protein are active 
during the cellular division (mitosis and meiosis, see e.g., H. Ris, 
1947). In cases of somatic mutation differences in both these sub- 
stances have been also encountered (cf. Caspersson and Santesson, 
1942). It has been suggested however that a major part of the amino 
acids of the tobacco mosaic virus is not responsible for its activity 
and for its reproductive ability (Schramm and Miiller, 1940; Miller 
and Stanley, 1941), whereas the nucleic acid seems to be essential for 
this purpose (Schramm, 1941). 

The theory here applied implies that four ionizations in any place 
within the: volume v produce the mutation of the virus. Since this 
volume (1.42 or 0.55 X 10° uw’) is quite large on atomic scale the 
question arises how can four ionizations anywhere within this large 
group of atoms produce such a fine effect as a mutation from one type 
of virus to another one. A possible explanation of this could be seen 
in the fact that the protein molecule contains many identical amino 
acid molecules so that a change in any of many of those molecules 
could produce the same effect. Another explanation which has been 
suggested in similar circumstances is that the energy has a great 
facility of transfer from one place to another within the volume 
(Timoféeff-Ressovsky, 1937, p. 1382). Quite convincing physicochemi- 
cal examples of ‘actual transfer of energy over volumes even larger 
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than this actually exist (Timoféeff-Ressovsky and Zimmer, 1947, pp. 
107-117; Scheibe, 1941a). More specifically with reference to the to- 
bacco mosaic viruses it is known that forces over distances of 100 A 
and more act between them (Stanley, 1939) ; the fact that these vi- 
ruses arrange themselves in a regular fashion even in weak solutions, 
so as to fill uniformly the available space with an interparticle dis- 
tance up to 500 A (Bernal and Fankuchen, 1941), is an excellent proof 
of forces acting between them over such distances, and where forces 
act energy is transmissible. Taking as the size of (v) even the larger 
figure 1.42 X 10° u?, each of the four ionizations would have to be 
effective within 0.355 X 10° A’, that is the energy available at each 
such event would have to travel in average perhaps not more than one- 
half of the diagonal of a cube of 0.355 X 10° A? which is only 35.4 
A. However because of the oblong shape of the protein and of the 
nucleic acid molecules of the virus the energy may have to travel over 
distances more than this. 

If instead of the whole virus, the nucleic acid alone were the site 
of the primary mutational process, a picturesque mechanism of 
transfer of energy would be possible. The determinations of the 
size and shape of the molecule of the nucleic acid of the tobacco mo- 
saic virus have shown that this molecule (700 A X 425 A?) is too long 
by a factor of several units to fit within the basic cell of the virus; 
nor can the nucleic acid molecules be arranged longitudinally along 
the virus, because their total length exceeds the available one by a 
factor of about 2 (Cohen and Stanley, 1942). The only possibility 
left is that the nucleic acid molecules are situated in a protruding 
fashion along the protein. This is confirmed by the fact that if the 
nucleic acid lie within the protein it would be recognizable there by 
means of X-rays because the absorption coefficient per unit of volume 
is more than two and one-half times as high for the nucleic acid as it 
is for the protein. Neither X-rays nor the electron microscope reveal 
the presence of such heavier molecules of the size of the nucleic acid 
within the virus (Cohen and Stanley, 1942). Further confirmation 
of this location of the nucleic acid may be seen also in the fact that 
an elimination of the nucleic acid from the virus by biochemical means 
(Schramm, 1941) does not split up its protein molecule. Now, if one 
imagines the nucleic acid molecules arranged in a parallel fashion 
with respect to each other (cf. Timoféeff-Ressovsky and Zimmer, 1947, 
p. 115; Butenandt et al, 1942), the virus would resemble a longitudi- 
nal bar to which several transversal bars are attached laterally ; since 
these bars would be elastic an impact on one of them would put the 
others in vibration. This analogy could visualize the transfer of energy 
along the virus when an ionizing particle collides with one of the pro- 
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truding nucleic acid molecules. The fact that the number of these 
molecules in a virus is 8 (Cohen and Stanley, 1942) and that the num- 
ber n of absorptions of energy sufficient for the mutation is 4 gives 
two possible pictures of the mechanism of mutation of the virus. 


FIGURE 2. Schematical diagram of a protein molecule with 8 transversally 
arranged nucleic acid molecules making up the tobacco mosaic virus. The dia- 
gram illustrates how 4 impinging ionizing particles (represented by circles) can 
uniformly affect all 8 molecules through energy transfers (arrows) over dis- 
stances not exceeding that of two closest neighbors. The theory implies that 4 
ionizations suffice for the induction of mutation. 


(a) If ionizing energy reaches a certain part of a nucleic acid 
molecule, a part of this energy is used for a structural change in this 
molecule and the remaining part is transferred to the neighboring 
nucleic acid molecules and possibly also to the protein molecule (Fig. 
2). Because of the long distance between the neighbors, some 600 or 
700 A, the transfer of energy would be limited to the closest neigh- 
bors, and the energy reaching a nucleic acid molecule by such trans- 
fer would not be sufficient to produce in it the entire change neces- 
sary for the mutation; however if a sufficient amount of energy is 
transferred also from the other neighbor the whole change of that 
molecule necessary for the mutation would be accomplished. Although 
this transfer would occur over quite a substantial distance of 600 or 
700 A, it should be kept in mind that transfers of energy over much 
longer distances are also known to occur in physicochemistry (Riehl, 
1941, p. 56; Scheibe, 1941b; Scheibe and Kandler, 1938; Méglich and 
Schon, 1938). This mechanism would imply that the mutation of the 
tobacco mosaic virus to the aucuba type requires a change in all the 
nucleic acid molecules of the virus. 

(b) If a change in one-half of the number of molecules of the 
nucleic acid were sufficient for the particular mutation considered no 
transfer of energy between these molecules needs to be invoked. 

The analysis can be carried out further by asking what part of 
the nucleic acid, if any, may be responsible for the mutation. If phos- 
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phorus is taken under consideration we obtain as the size of the 
volume (v): 


vV=5.5 X 107 X 12.46 X 1.6/150 = 7.31 X 10-8 ae, 


because 150 cm~ is the volume absorption coefficient of phosphorus 
at 1.5 A. But the volume of phosphorus existing in the whole nucleic 
acid of the virus is 


2.23 X 10° X 1.6 X 0.0913/2.2 = 14.8 X 10 13, 


which is twice as much as v. Consequently, changes in one-half of 
the phosphorus existing in the nucleic acid of the tobacco mosaic virus 
would be sufficient to produce its mutation to the aucuba type. Such 
changes could involve for instance the phosphorus of every second 
molecule of the nucleic acid and would not require transfer of energy 
between various molecules of this acid. According to J. Johnson 
(1949) there are some twenty variants of the tobacco mosaic virus. 
Theoretically one cannot exclude the possibility that mutations to 
other variants imply the same type of mechanism, different molecules 
of the nucleic acid being involved in each case. 

The table below summarizes the alternative conclusions of the 
analysis. 


MUTATION OF THE TOBACCO MOSAIC VIRUS TO THE 
AUCUBA TYPE INDUCED BY 1.5 A X-RAYS. 
POSSIBLE SITES OF THE PRIMARY MUTATIONAL PROCESS 


Phosphorus of the 
Nucleic Acid Nucleic Acid 


5.5 X 10-7 Gal Se als 
Te 8 


Site (S) Protein 


Size of the 1.42 X 10° 
volume (v) us 
within which 

the primary 

process occurs 


Approximate 35 4 
ratio of the 

volume of (S) 

to the volume 

of (v) 


(S) is a basic 
hexagonal cell 
unit as deter- 
mined by X-Ray 
diffraction 
(Bernal and 
Fankuchen, 1941) 


Possible in- 
terpretation 


At least two 
molecules of the 
nucleic acid are 
affected by the 
process 


Phosphorus in 
at least every 
second molecule 
of the nucleic 
acid is involved 
in the mutation 


N. B. Calculations rade on the basis of 5 X 10-5 »3 as the volume of the virus. 
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As much as this theory seems to agree with all the pertinent 
known facts about the tobacco mosaic virus there may be perhaps 
some doubt whether neglecting the inactivation of the aucuba viruses 
by the X-rays after they have been formed by this same radiation 
from the tobacco mosaic viruses is not a too crude approximation. In 
fact an inspection of the data of J. W. Gowen (1941) shows that the 
number of aucuba lesions obtained from an irradiated sample of the 
tobacco mosaic virus increases with the dosage up to a certain limit 
only to decrease thereafter. A possible explanation of this fact is the 
inactivation and to a small extent perhaps also mutation of some of 
the aucuba viruses which were previously formed from the tobacco 
mosaic viruses in the course of the same experiment. These second 
order effects make the number A in equation (7) appear smaller than 
it would have been if one could take them into account. We will see 
later that these effects are counterbalanced in part at least by the 
concentration effect. However before we consider the latter we pro- 
ceed to a type of analysis in which different kind of second order 
effect makes A appear larger than it actually is. 

(ii) Inthe previous analysis the inactivation of the aucuba vi- 
ruses formed during the experimentation was neglected. An oppo- 
site type of approach is obtained by assuming that equal fractions of 
both the original tobacco mosaic viruses and the newly formed aucuba 
viruses are inactivated during the experimentation. Such approxima- 
tion will make A/T appear larger than the value this ratio should have 
to correctly express the efficiency of X-rays as the mutating agent of 
the tobacco mosaic virus to the aucuba type. In fact the mutated vi- 
ruses in their new aucuba form are obviously subject to a smaller 
amount of radiation than the original tobacco mosaic viruses. The 
implications of this assumption will now be analyzed. As to the re- 
verse and other possible mutations of the aucuba viruses, we will 
neglect them as in (i); this can not change our conclusion that A/T 
appears here larger than it should because mutations are induced at 
a much smaller rate than inactivations. 

The analysis is carried out here in a similar way as in (i) with 
the only difference that no survival probability needs to be consid- 
ered since this is now assumed to be the same for both tobacco mosaic 
and aucuba viruses; instead of (7) the following relation is used now: 


P(Ng + Np) =A/T. 


Figure 1 and the table give the results of calculations for n = 1 baad 
3. The standard deviation increases with n. An inspection of the 
figure shows also that the fit cannot be improved by increasing the 
value of n , because the experimental data indicate a pronounced con- 
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cavity downward whereas the theoretical curves have the opposite 
concavity for n > 1. In fact, formulae (9)—(10) show that the re- 
lationship between P and N, is the incomplete gamma function 
y (UN. D n), except for the scale of ordinates and for a displacement of 
the origin of the system of coordinates. Because of the smallness of 
P(Ns), such displacement does not change the concavity of the curve 
P(N>) in the present range at the origin. [For a graphical represen- 
tation of the curves y see e.g. P. Curie (1929), J. H. Webb (1948), J. 
Th. Van der Werff (1948), p. 135.] It is seen that the best fit that 
can be obtained with this theory is forn = 1. It is not satisfactory 
however since it gives a standard deviation for v of 15.7%, conse- 
quently on the basis of the latter the choice between the two analyses 
is easy since analysis (i) implies o = 2%. There is also another very 
important reason in favor of analysis (i): To produce the lesions the 
plant is infected by spreading the virus on the leaves. Now when the 
virus is used in high concentration there will be many cells which 
more viruses enter than it is necessary to produce a lesion (Caldwell, 
1933). In the experiments of J. W. Gowen (1941), which are here 
analyzed, the tobacco mosaic virus was used in high concentration, 
whereas the aucuba virus obtained through mutation was in a very 
low concentration since most of the viruses of this sample were inac- 
tivated by the X-rays; consequently, as far as the concentration effect 
is concerned, A in equation (7) has its correct value corresponding 
to the number of aucuba viruses existing in the sample, whereas T 
appears smaller than it should. In this way A/T of that equation ap- 
pears larger than it would have been if the concentration effects 
have been taken into account. But we know from a previous discus- 
sion of analysis (i) that 4/T appears also smaller than it should be- 
cause mutations and inactivation of aucuba viruses have been neglect- 
ed in that analysis. In conclusion we may say that the approxima- 
tions used in analysis (i) counteract and in analysis (ii) accentuate 
each other. We see in this way that we have all the reasons to accept 
the results of analysis (i) against those of analysis (ii). 

Both types of analyses (i) and (ii) make it possible to calculate 
Ng from equation (11) after v and have been determined. Using 
the results of analysis (i) we obtain Nx = 282000 or 729000 ioniza- 
tions per 2 according to whether (v) is situated in the protein or in 
the nucleic acid. If the result of analysis (ii) were accepted it would 
give Ns = 72400 wr. Two alternative interpretations of these figures 
are possible. Let us consider for instance Ns = 500000 as the average 
of the two figures of analysis (i). The cosmic radiations and the 
y-rays from rocks amount to about 0.2 roentgen per year (see e.g. 
Compton, 1942), this gives in a round figure 1/2 ionizations per u 
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per year. Consequently if the spontaneous mutation were the exclu- 
sive work of naturally occurring ionizing energy the nature would 
have to work one million of years to do it! On the other side if the 
spontaneous mutation were carried out entirely in a single generation 
of the virus, mutating agents should naturally occur in the life of the 
virus more than one million times as powerful as the naturally existing 
ionizing energy, since the length of a generation of a virus is a very 
small fraction of one year. It is difficult to see which one of these in- 
terpretations is closer to the truth, since our understanding of spon- 
taneous mutations is very incomplete. 

In the numerical calculations outlined in analyses (i) and (ii), the 
tables of the incomplete gamma function of E. C. Molina (1942) have 
been used. The tables of K. Pearson (1922) are less convenient be- 
cause of a transformation of the argument which they involve. Since 
the mathematical problem in (i) and (ii) consists of calculating ¢ 
from assigned values of n and of y= y(c, n), the expansion of G. A. 
Campbell (1923) for the inversion of the incomplete gamma function 
by means of the probability integral could be also used instead of the 
tables of the y-function. 

The two analyses presented here each imply an opposite type of 
approximation; in the first analysis the inactivation of aucuba viruses 
formed during the experiment was neglected; in the second analysis 
the percentage of the same aucuba viruses inactivated has been as- 
sumed to be equal to the percentage inactivation of the original to- 
bacco mosaic viruses. The exact situation lies somewhere between 
these two approximations. Which of them is better, it is difficult to 
say a priori. The first analysis shows a very outstanding agreement 
with the known experimental facts; this may be due to the concen- 
tration effect which counterbalances the approximation used in this 
analysis. 

A mathematically more detailed type of analysis is outlined in 
the appendix. The paucity of experimental data does not justify car- 
rying out this analysis more extensively. 


Appendix 

An analysis which takes into account the different average dos- 
ages which the two types of viruses must survive in order that the 
mutation be detected could be the following. Let exp(— k x) be the 
probability of a tobacco mosaic virus surviving as such a dose x, and 
exp (— h x) a similar probability for the aucuba virus. The expo- 
nential form of these probability functions is suggested to some ex- 
tent by experimental data [cf. formula (6)]. If a tobacco mosaic 
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virus has survived as such a dose x its probability of mutation to au- 
cuba with an additional dose dz is according to equation (5) 


d=e" (vx)™ 1d (av) /(n—1)! 


Consequently the probability of mutating a tobacco mosaic virus with 
a dose < Np and of detecting such mutation with the total dose = Np 
is 


PN 6) i d-exp [— (Np—2z)h] 


a0 

12 

=vu(vt+tk—h)"y[(v+k—h)Ny,n] exp(—Nph). or 
If the data on inactivation of the aucuba and the tobacco mosaic 

viruses of J. W. Gowen (1940, 1941) are plotted against the dosage 

quite a substantial dispersion appears. The expression 


S=exp(— 5 X 10-* D) Gis) 


for the probability of survival of each variant is about the best of 
exponential type. The dose D is here expressed in roentgen. Since 
this relation implies h =k, equation (12) becomes 


P(Np)/S(No) = y (UN, Nn). 


Taking into account the spontaneous mutation and assuming that 
there is no spontaneous inactivation the previous relation should be 
written: 


P(Np + Ns) /S(No) = y(UNp + UNsg,). 


Now, if the left-hand side of this equation is plotted against D using 
for S formula (13) and for P(N, + Ns) the ratio Ap/T, where T> 
is the number of the tobacco mosaic lesions obtained before adminis- 
tration of X-rays and A, is the number of aucuba lesions obtained 
after administering a dose D to the tobacco mosaic virus, the curve 
obtained shows a pronounced maximum. This is incompatible with 
the form of y which is a monotonously increasing function. The conclu- 
sion is that the assumption h =k, in particular relation (13), is a 
too crude approximation. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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After some general remarks the results of the approximation method 
of Kryloff and Bogoliuboff are applied to show the existence of a thresh- 
old for a simple non-linear differential equation of a special form. 


Many oscillation phenomena can be described, at least as a rather 
good approximation, by the differential equation 


ttaz+p=0 (1) 


in which the time-dependent variable « represents the deviation from 
the equilibrium position and a and f are constant with respect to the 
time ¢ as well as to x. If a is positive the corresponding oscillation 
is damped and it dies out, while when a is negative the damping is 
negative and the amplitude of the deviation increases with the time ¢. 

When the coefficients a and 6 are themselves functions of the 
variable x , the problem is somewhat more difficult. B. van der Poll 
(1930) has discussed the solution of the equation 


a—e(1—2)e+a=0. (2) 

In the case ¢« << 1 he has shown analytically that for large values of 
t the solution of (2) becomes 

x=2sin(t+ 4), (3) 


quite independent of the initial conditions; it represents a stationary 
or self-excited oscillation (van der Poll, 1930, Fig. 10). For the case 
« > 1 van der Poll used the method of isoclines; he replaced (2) by 
the set of differential equations 


£=y; 
(4) 
y=e(1—2*)y— 2; 
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_, ay 
and plotted in an (x,y) diagram the curves along which de = con- 


stant. By this method he noticed that the integral curves of equation 
(2) approach a fixed closed curve, quite independent of the initial 
conditions. These “limit cycles” of Poincaré represent also in this 
case a stationary periodic solution. 


The coefficient of x in (2) is not constant, but a function of 2, 
which is negative for values of x between + 1 and —1, so that for 
those values of « the damping is negative. One might expect from 
this that the differential equation 


a+e(x—p)(x—qxa+u=0, (6) 


in which p and q are both positive constants such that p < q, would 
also give a stationary periodic solution, independent of the initial con- 
ditions. If this was the case, the simple equation (6) would describe 
a threshold phenomenon. To illustrate this, notice that the coefficient 


of x in (6), and therefore the damping of the motion described by 
(6), is positive for « < », but negative for p < x < q and again posi- 
tive for x > q. Therefore we get in an (2, ¢) diagram the next re- 
gions as shown in Figure 1. If the initial conditions are such that x 


xX 
REGION OF POSITIVE DAMPING III 
q SS 5 SS a SS SS Se 
REGION OF NEGATIVE DAMPING II 
p ae 


REGION OF POSITIVE DAMPING I 


FIGURE 1 


remains smaller than p, the oscillation is always positively damped 
and dies out. But if the initial conditions are such that the point in 
Figure 1 which represents the state of the system comes in region II, 
there might be a possibility, just as above with the negative damping 
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in the case of van der Poll, that a stationary oscillation, quite inde- 
pendent of the initial conditions, thus self-excited, exists. By com- 
bining the two last mentioned cases, we notice the possibility of the 
existence of an all-or-none phenomenon in the case of equation (6). 
However this is not the case. Though the only difference between (6) 
and van der Poll’s equation (2) is that » in the former and in Figure 
1 has changed its sign, this gives an essential difference. To show 
this we apply the result of a method developed by J. Sjohat (1944). 
He considers the more general equation: 


a—eF(x)n+nr=0. (7) 


Assuming that this equation has a periodic solution «(t) with period 
T, he expands this function x(t) in a Fourier series. Because 


fe 
af cdt=0, 


as can be seen by integrating (7) from t = 0 to t = T, the Fourier 
expansion of x(t) can be written in the form 


x(t) =D A, sin (nyt +1d,) , (8) 
Sjohat shows that the first term in this expansion gives a good ap- 
proximation of x(t). Therefore 


x(t) =A, sin(vyt +44). (9) 


Moreover he evaluates A, for the case in which F(x) is a poly- 
nomial: 
F (a) =65 Fiat ess oe’ « 
Under the assumption that equation (7) has a periodic solution « (¢) 
and for this form of F(x) he finds that z = A,? is approximately a root 
of the next equation: 


See Rape oe ek, (10) 
Bee a sabi Ob | 


ge? 
in which m is the largest integer ins 4 


From this equation we see that equation (6) has no stationary 
oscillations as solutions, because in the case of that equation r= 2, 
m = 1 and a and a have the same sign. Therefore the necessary con- 
dition that (10) has at least one non-negative root cannot be satisfied. 
The reason for this failure becomes apparent. by comparing this case 
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with that of van der Poll. In the latter case there is also a region of 
negative damping for negative values of ~, namely that between 0 
and — 1, while in the former case the damping is positive for all 
negative values of x. Therefore the movement, when the deviation 
x has become smaller than p and afterward, becomes negative; it is 


€=01 


-20 


-05 
-1.0 


P - Ay 4) —20 


2.0 


-2.0 


05 


Ficure 2. Integral curve of the equation 
& + e(%—1)(x—3)z+"—=0 


having the initial conditions x, = —4 and x, = 4 for the case «= 0.1. 
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so strongly damped that x does not exceed p again, even when it has 
done so once before. We have confirmed this presumption by using 
the method of isoclines (Fig. 2). 


This explanation of the negative result has lead us to expect 
to get the description of a threshold phenomenon if the coefficient a 


of x in equation (1) is also negative for some negative values of x in 
such a way that for negative values of x also we have a situation as 
in Figure 1. In other words, we extend Figure 1 to the next figure. 
(Fig. 3). We will show that this is the case by examining the example 


X 
REGION OF POSITIVE DAMPING IIT 


REGION OF POSITIVE DAMPING III’ 


FIGURE 3 


in which « in equation (1) is the next function of x: 
a=e(4—1)(*#+1) (x—38) (4+ 8), 
or 
a= e(z*— 1027+ 9). 


Therefore, we will investigate the differential equation: 


a+ e(st—10a2+9)a+a=0. (11) 


The graph of this function a of x has the form given in Figure 4. We 
will apply the method of Kryloff and Bogoliuboff (Minorsky, 1947). 
By this method it is shown that in the case of a more general equation 
than (6), namely 
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FIGURE 4 


etef(x,2) +ox=0, (12) 


in which case « << 1, the time-dependences of the amplitude a and 
the phase 6 = wt '+ ¢ of the approximate solution 


x=asin(wt+ ¢), 


are given by 


ee =-“= f° f(asin py, aw cos p)cos pdy, (13) 
dt aw 2n 


and 


oe f(asiny, aw cos y)sin ydy. 
dt aw 2a vo 


By applying (13) to our case of equation (11) in which we have 
o—=1and 


i dl ree “) = (a*— 102? +9)2, 
therefore 
f(asin yp, aw cos yp) = (a* sint y — 10a? sin? y + 9)aw cosy. 
We then get 


Cae tele aes ‘ 
— = — — (a4 sin‘ y — 10a? sin? y + 9) aw cos? pdy, 
dt w 22 Jo 
or (15) 
da 1 5 9 
So Se (=ga Fore) 
dt 16 4 2, 
For a steady state we have 
da _ 
the ae 


and therefore in this case 
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: a> ae ale 0 
ae — 6 te 0, 

16 4 2 
Hence a,—0 » 2,3 = + 2.17, a5 = + 3.91. Therefore there are two 
stationary oscillations possible, one with amplitude 2.17 and the other 
with amplitude 3.91. 


From the following form, in which we can now write (15), 


ot (gt — 2.17) (# — 3.91%) 
dt 16 ‘ ; 
we see that 
da 
Por 0 <¢ 2 217 —<0 
dt 
da 
KonZaliqea—3.08 Fee 
da 
For a > 3.91 —<0. 
dt 


Therefore there is in the case of equation (11) a stable stationary os- 
cillation, which has the amplitude 3.91. 
The frequency 2 of that oscillation can be calculated from (1) 
according to 
PECbts de 


Oo; 


Q=—= : 
dt dt 


In the case of equation (11) this gives: 
1 2 : 
Q=ar+ ~—f{ (a* sint yp — 10a? sin? y + 9) aw cos yp: sin yp dy 
aw 2n J, 


=o=1, 


the integral being zero and w = 1 in equation (11). Therefore the fre- 
quency of the oscillation in the case of equation (11) is constant, as 
for the van der Poll equation (2). However in general the frequency 
will depend on the amplitude of the oscillation. 

I wish to thank Professor N. Rashevsky for the suggestion of 
this problem and his generous assistance. 

This work was aided by a grant from the Dr. Wallace C. and Clara 
A. Abbott Memorial Fund of The University of Chicago. 
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The fundamental idea in the background of this book is derived from some 
purely formal quantitative similarities between phenomena of electrical nerve ex- 
citation and radiobiological effects. Because of the complexity of these phenomena 
and our very rudimental knowledge of their biophysical nature it is difficult to 
say how much basic truth underlies such similarities. It is certain that the ionized 
state of matter is the primary intermediary of both processes. But to what extent 
is this common factor of importance? Certainly the energies involved are in both 
cases of a completely different order of magnitude; so far no nerve excitation by 
ionizing radiation has been reported, and in the clinical field electrotherapy and 
radiotherapy reveal completely different physiological mechanisms. The thresh- 
old concept, so firmly established in the electrophysiology of the nerve, seems to 
have a somewhat less fundamental importance in radiobiology. In fact, whereas 
for instance the erythema threshold is an accepted fact, the induction of mutations 
in Drosophila by X-rays, a most carefully investigated field of radiogenetics, does 
not reveal any threshold phenomenon. 

R. M. Sievert (1941) was the first to formulate mathematically the existence 
of a reserve of a certain substance which has to undergo a certain change before 
the radiation may manifest its biological effect. If the energy administered is in- 
sufficient to carry out such a change no effect would appear. Along the same line 
of thought are the current suggestions that some biological effects of radiations 
consist of an inactivation of the existing supply of enzymes. The amount of this 
supply would then determine the magnitude of the threshold, if it exists. 

The author starts his book by pointing out some formal analogies between: 

(a) the Bunsen-Roscoe or Schwarzschild law and the relationship between 
the intensity of a constant current and its minimum time of application to cause 
a muscle contraction; 

(b) the latency of several radiobiological phenomena and the persistency of 
certain effects in the nerve after the cessation of the cause as evidenced, for in- 
stance, by a succession of subliminal stimuli inducing nerve excitation ; 

(c) the oscillatory aspects of many phenomena of electrical nerve excitation, 
for instance the change in excitability due to short subliminal condenser dis- 
charges, and the oscillatory character of certain radiobiological processes stressed 
so much by A. Forssberg (1941, 1943) in his remarkable experiments on the ef- 
fect of extremely small doses of X-rays (one thousandth of one roentgen) on the 
growth rate of the fungus Phycomyces Blakesleeanus. 

This oscillatory character is known to occur in many other effects of radia- 
tion, of which the common erythema wave is perhaps the oldest although not the 
most evident example. Reporting these oscillations in the growth rate of irradi- 
ated root tips of Vicia faba, L. H. Gray (1949) suggested that this effect may be 
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apparent only due to some geometric peculiarities of growth of the irradiated 
material. However oscillatory character of radiobiological phenomena is known 
to exist in many other cases, and has been observed also in purely radiochemical 
processes such as irradiation of colloids (J. A. Crowther, H. Liebmann, R. Jones, 
C. C. Mills, 1936-1940). The recent experiments of B. E. Proctor and S. A. Gold- 
blith (1949) on the change in the ultraviolet absorption spectrum of apricot ex- 
tract and of 5- (hydroxymethyl) furfural with doses up to 20 million roentgen of 
cathode rays also show a clear indication of a beginning oscillation. The mechan- 
ism of such oscillations in radiobiology is not known; although suggestions for its 
explanation have been made (eg. A. Forssberg). Oscillations in diffusion with 
intervening chemical reactions have been studied mathematically by several au- 
thors (D. G. Kendall, 1948; N. Rashevsky, 1948) and an attempt of their applica- 
tion to those still non-understood oscillatory radiobiological phenomena would 
seem to be a worthwhile enterprise. 

Although the previously published work by the author lies in the field of 
theoretical radiobiology, two-thirds of the present book is neurophysiology. It is 
in the above-mentioned experiments of A. Forssberg and in the mathematical work 
of R. M. Sievert (1941) that the book had its origin. Both Sievert and Forssberg 
had a clear notion of the latency of the effects observed by them; this, combined 
with the assumption of an active reserve or ability of the organism to recover lead 
them to an interpretation of the oscillatory character of their experimental re- 
sults. Following essentially the line of thought of these Swedish radiologists, the 
author assumes that the phenomena in both fields, the neurophysiological and the 
radiobiological, depend on the amount X of a certain substance ruled by a differ- 
ential equation of the type 


dX (t)/dt =—I(t—a) X(t) + L(t—}6)Y(‘), (1) 


where Y is understood essentially as the amount of the reserve substance, and J 
and L are coefficients describing the kinetics of two time-lag processes, as the form 
of the functions I(t — a) and L(t — b) indicates, t being the time and a, 6 the 
time-lags. 

The neurophysiological part of the book is worked out mainly under the con- 
dition Y = constant, which is Sievert’s assumption. The two possibilities a > b 
and b > @ are separately discussed. The author assumes 


I(t) = (14 ¢,i)I,; L(t) = (1+41)L,; (2) 
at the cathode, and 


Lee Lt ies eee, 
Li-- 0,2 1. ¢,4 


at the anode, where 7 is the intensity of the current and the c;’s, I, and L, are 
constant. The author looks upon I(t) and L(t) as expressing the concentration of 
cations. He defines L(t)/I(t) as the physiological electrotonus of the nerve. This 
definition plays an important part in the author’s systematics of his theory; a 
large part of his study concerns the case in which L(t)/I(t) is independent of i. 
Taking the smaller of the two quantities a, b as zero, which is equivalent to meas- 
uring the time from a suitable moment, the author’s equation reduces to either 
one of the following two forms: 


(3) 


dx (t) /dt = — I(t —a) x(t) '+ I(t); (4) 
da (t)/dt =—I(t)a(t) + I1(t—a). 
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The classical topics of electrical nerve excitation are discussed on the basis of 
these two equations. Since the latter are mixed differential and finite difference 
equations, they split up into four differential equations in the case of a constant 
current of finite duration. The functions x(t) represent then a wave-like curve 
with three discontinuities of its tangent line. If the current lasts indefinitely a 
similar curve is obtained with «—0 only at t= © and one discontinuity of the 
tangent line. The relationship between the minimum duration of a constant cur- 
rent and the threshold intensity of the latter comes out to be the same as in Hill’s 
theory. The author passes then to discuss the problem of condenser discharge 
through the nerve: he obtains essentially the same type of curve as with a con- 
stant current of finite duration, except that only one discontinuity of the tangent 
line appears. Here again the author obtains Hill’s relation between the threshold 
voltage and the capacity and resistance. He then applies his equations to the ac- 
commodation phenomena in the case of linearly and exponentially increasing cur- 
rents and discusses the problem of excitation at the anode by interrupting a 
constant current, under an arrangement in which a make excitation at the cathode 
is prevented. The relationship between the minimum intensity of the current and 
its duration necessary to obtain a break-excitation gives a finite duration for an 
infinite current which does not agree with Hill’s formula, although the author 
contends that the experiments of Cardot and Laugier (1912, 1913) suggest the 
correctness of his mathematical result. He discusses then the case of a constant 
current with a gap. In the case of alternating currents the threshold intensity 
comes out directly proportional to the frequency for high frequencies and inverse- 
ly proportional for low frequencies. The author stresses the possibility of de- 
scribing by his equations the experiments of J. Erlanger and E. A. Blair (1981) 
on the change of excitability after a subliminal direct current has been put 
through the nerve, without discussing the question in detail, however. The ex- 
periments of W. A. H. Rushton and B. Katz on the decline of nerve excitability 
after short current pulses show some deviation from the author’s theory. He dis- 
cusses also A. Chweitzer’s measurements (1937) on the change of rheobase 
after the passage of a current of constant intensity. His calculations do not give 
a minimum as flat as that obtained experimentally, however the probabilistic ap- 
proach which he considers later leads to a better fit. , 
The author passes then to a brief treatment of the case in which his physio- 
logical electrotonus depends on the intensity of the current, that is when the ¢;’s 
in equations (2) and (3) are different from each other. This approach leads in 
many cases to no essentially new conclusions. Since a careful discussion of the 
two cases a S b gives fundamentally the same results, the author sees no way of 
deciding which one of the two inequalities actually holds. He seems to be particu- 
larly concerned with the following contradiction between experimental findings 
and a consequence of the theory presented by him in the first hundred pages of the 
book: the ratio \ of the rheobasic current i to the minimum value of di/dt causing 
nerve excitation comes out to be smaller than a or 6 in the case of an exponential- 
ly increasing current, whereas from an analysis of the experimental data it fol- 
it should be larger. ; 
Se seenree this and ee discrepancies the author generalizes his eer 2 
assuming that the functions I(t) and L(t) which characterize the news ile o 
bility represent cumulative effects of a number of random events accor ee 0) 
Poisson distribution. Under this assumption the mathematics becomes quite com- 
plicated. Some of the relations which the author obtains by this new theory, ine 
far as numerically analyzed by him, do not differ qualitatively from the preceding 
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ones, except that no angular points appear on his curves. Quantitatively, how- 
ever, there are cases which give by this theory essentially different results and 
this seems to be the reason why the author considered this modification of his pre- 
vious approach. For instance, prevention of excitation after an overthreshold cur- 
rent pulse has been applied, by means of a strong current pulse in the opposite 
direction immediately after the first pulse (Rushton and Katz) can be described 
by this theory but not by the preceding ones. Similarly Hill’s factor \ is now not 
necessarily smaller than the rheobasic utilization time, as the previous theory in- 
dicated. The theory of Chweitzer’s experiments now gives a flatter minimum than 
before. The author claims that this is the only mathematical account of Chweit- 
zer’s results known up to now. 


The last part of the book is entitled “Biological Reactions Caused by X—and 
y—Rays.” It is not clear to the reviewer why the author restricts the title to these 
two forms of radiation only. Since his theory is descriptive and formal there 
seems to be no reason for discrimination between the various types of radiations, 
including ultraviolet and even the visible range. The author’s paper of 1942 was 
an expansion of Sievert’s mathematical concepts and the radiobiological part of 
the present book is a further development of these ideas. A whole chapter deals 
with an interpretation of Forssberg’s experiments by means of an equation of 
type (1). The remaining chapters discuss other radiobiological experiments re- 
vealing oscillatory character. In considering the relation 7t? = constant hetween 
radiation intensity 7 and the time t of its application necessary to achieve a certain 
effect, the author interprets the case p <1 as due to recovery phenomena which 
are more efficient at lower radiation intensities if the total administered energy is 
the same. He interprets the case p > 1 in the following way: 


With an irradiation of very short duration only certain cells, being in a 
very radio-sensitive phase, may be killed. With an irradiation of longer 
duration there may be more chance that other cells come into that radio- 
sensitive phase and so the effect may be enlarged by diminishing the 
intensity. 


If the radiation acts through a catalyst, the effects of radiation will depend not 
only on the amount of this catalyst but also on the time during which it is active. 
This interpretation of p > 1 has been suggested by the reviewer (1947) as a pos- 
sible explanation of the relation it? = constant found by A. M. Brues, H. Lisco 
and M. Finkel (1946) for the induction of bone sarcoma in mice by a radioactive 
isotope of strontium. It should be kept in mind however that no direct evidence 
of the existence of such catalyst has been obtained up to now. 


The book closes with brief discussions of the Bunsen-Roscoe law, the cumu- 
lative dosage concept, the F. Dessauer—M. Blau—K. Altenburger approach to 
radiation effects through the cumulative Poisson distribution (1922). As a whole 
the radiobiological part appears to be treated in a rather sketchy form, if one 
compares it with the systematic and quite detailed mathematical work presented 
by the author in the part dealing with the electrical nerve excitation. The reason 
for this lies probably in the much more pronounced complexity of radiobiological 
phenomena as compared with the relatively simpler situations met in exciting: 
electrically a nerve in vitro. Contrary to the author’s impression he uses mathe- 
matics for a quantitative description and not for a physical explanation of the 
observations. A comparison of the radiobiological part of his book with D. E. 
Lea’s “Actions of Radiations on Living Cells” may give an idea to what extent:the 
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reviewer’s viewpoint is correct. However the history of science shows that a new 
theory starts sometime with formal, somewhat abstract and purely descriptive 
concepts before acquiring a concrete reality in terms of experimentally known 
facts. From this viewpoint the book may be a very interesting one, and is certainly 
suggestive of many ideas deserving further thought. 


I. OPATOWSKI 
The University of Chicago 
Committee on Mathematical Biology 
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